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SECTION  I 


INTRODUCTION 

This  report  describes  the  results  of  theoretical  and  experimental 
investigations  of  flame  temperature  methods  based  upon  laser  excited 
fluorescence  measurements  of  inorganic  probes  introduced  into  the  flame 
gases;  these  studies  were  conducted  under  USAF  Contract  F3361S-78-C-2038 

The  objective  of  this  contract  was  to  investigate  the  use  of  laser 
excited  fluorescence  for  spatial/temporal  temperature  measurements. 
Development  of  such  new  combustion-diagnostic  techniques  is  important 
to  the  Air  Force  in  connection  with  the  development  and  refinement 
of  cleaner  and  more  efficient  jet  engines. 

Atomic  fluorescence  excited  with  a  spectrally  quasi-continuum 
source  has  been  shown  to  be  a  reliable  method  of  measuring  flame  tem¬ 
peratures.^  In  fact,  it  combines  the  advantages  of  optical  methods 

i.e.,  of  not  disturbing  the  combustion  process,  with  the  capability 
of  providing  spatially-resolved  measurements,  the  last  feature  being 
common  to  all  "scattering  based  methods",  such  as  for  example,  the 

Q 

Raman  method. 

The  use  of  the  laser  as  an  excitation  source  for  atomic  fluores¬ 
cence  has  the  obvious  advantage  that  spatial  resolution  can  be  signi¬ 
ficantly  improved  because  of  the  low  divergence  of  the  beam  coupled 
with  its  high  spectral  irradiance.  The  use  of  a  pulsed  laser/gated 
detector  system  also  has  the  advantage  of  temporal  resolution. 


Moreover,  as  shown  by  Measures^  and  Omenetto  and  Winefordner ^ 
it'  the  laser  irradiance  is  such  that  saturation  of  the  excitation 
transition(s)  can  be  achieved,  then  the  temperature  of  the  atomic  sys¬ 
tem  can  be  evaluated  without  the  necessity  of  calibration  of  the  electro- 
optical  detection  system  and  without  measuring  the  laser  power.  How¬ 
ever,  as  discussed  below,  other  constraints  are  imposed  on  the  method. 

The  aim  of  this  report  is  to  give  for  the  first  time  several  means  (5 
different  approaches)  of  calculating  the  temperatures  of  several 
based  laboratory  flames  by  laser  excitation  of  a  3-level  atomic  probe, 
such  as.  In,  Tl,  Pb ,  etc.,  aspirated  at  a  low  concentration  into  the 
flame  and  excited  simultaneously  or  sequentially  with  one  or  two  laser 
beams . 
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SECTION!  II 


THEORY  OF  FLAME  TEMPERATURE  METHODS 

The  theoretical  steady  state  expressions  have  been  derived  by 

1 2 

the  use  of  the  radiance  expressions  derived  by  Boutilier,  et  al 

1 2 

In  Appendix  A,  (the  Boutilier,  et  al  paper  ),  more  details  are 
given.  To  facilitate,  the  discussion,  it  will  be  assumed  that  the 
atomic  probe  introduced  into  the  flame  gases  has  3  significant  levels, 
where  level  1  is  the  ground  state,  level  2  is  the  closely- spaced 
metastable  state,  and  level  3  is  the  excited  state.  In  Figure  1,  this 
representation  is  given  along  with  the  significant  rate  constants; 
k's  are  radiationless  pseudo  first  order  rate  constants;  A's  are  radia- 
tional  rate  constants;  B's  are  Einstein  induced  absorption  or  induced 
emission  coefficients:  Ev ' s  are  source  spectral  irradiances;  and  c  is 
the  speed  of  light.  A  further  discussion  of  these  terms  is  given  in 
Appendix  A.  In  Figure  2,  the  specific  energy  level  diagrams  of  Ga, 

In,  and  T1  are  given.  In  the  following  textual  passages,  the  5  basic 
methods  for  flame  temperature  measurements  based  upon  laser  excited 
fluorescence  will  be  described  (also  see  Appendix  B) . 

Method  1.  Linear  Two  Line  Method  (based  on  method  in  Omenetto, 
et  al ^) . 

In  this  case,  the  ratio  of  fluorescence  radiance  signals  resulting 

when  exciting  1-3  and  measuring  the  fluorescence  of  3-2,  Bp  and 

r  3+  2 
1-3 

resulting  when  exciting  2-3  and  measuring  the  fluorescence  of  3-1, 

BP  .  The  resulting  temperature  expression  is 
^3-1 
2-3 


3 


Figure  1.  Energy  Level  Diagran  and  First  Order  Rate  Constants  for  Activation- 
Deactivation  of  a  Three  Level  Tl-Like  Atom. 
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where  is  the  excitation  spectral  irradiance  at  the  designated 
transition  (in  J  s'  nT^Hz”  )  ,  X  is  the  wavelength  (in  nm)  at  the 
designated  excitation  lines,  and  e12  1S  the  energy  of  the  metastable, 

2,  level  (in  eV)  .  This  method  of  flame  temperature  measurement  is 
independent  of  quenching  processes  but  requires  calibration  of  the 
spectrometr ic  system  to  enable  measurement  of  (relative)  values  of 
Ev  and  Bp  at  the  appropriate  wavelengths.  In  addition  to  this  approach, 
several  other  disadvantages  are  present:  (a)  changes  in  the  quantum 
efficiency  due  to  quenchers  in  the  flames  leads  to  a  change  in  Bp 

values  and  this  can  deteriorate  the  S/N  ratio;  (b)  the  source  spectral 
irradiance  must  remain  on  the  linear  part  of  the  Bp  vs  Ey  curve; 

(c)  at  low  E^  values  and/or  in  highly  quenching  flames,  Bp  will  also 
be  low  causing  S/N  problems;  (d)  the  post  filter  effect  leads  to 
deterioration  of  the  fluorescence  signal  especially  for  the  3-1 
transition.  Because  of  these  serious  difficulties  and  our  aim  to 
develop  methods  useful  for  combustors,  we  spent  most  of  the  contractual 
period  involved  in  the  development  of  methods  with  fewer  difficulties 
and  greater  chances  of  success.  Experimental  results  for  flame  tem¬ 
peratures  via  this  method  will  be  given. 


Method  2.  Saturation  Two  Line  Method  With  Sequential  Pumping 
In  this  case,  atomic  fluorescence  of  an  atomic  probe  is  produced 
at  3-*l  or  at  3-2  after  excitation  at  1-3  and  at  2-3,  respectively.  The 
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general  expression  for  the  flame  temperature,  T,  is 


E12/kT 

assuming  e”  <<  1,  and  letting  i  and  j  be  either  levels  1  or  2. 

The  terms  are  defined  as  follows:  Bp  is  the  same  as  before;  k's  are 
pseudo  first  order  rate  constants;  A's  are  Einstein  coefficients  of 
spontaneous  emission  (first  order  radiational  rate  constants),  u's 
are  frequencies  of  transitions,  and  g's  are  statistical  weights  (degen¬ 
eracies).  If  k21  >>  (A^p  +  kjp)  £  (A^  +  k32^  anc*  substituting 
for  the  statistical  weights,  then 


Finally,  if  n2  ^  30  n^  (n2  =  population  density  of  state  2  and  n^,  = 
total  population  density  of  all  levels),  then 


The  advantages  of  this  approach  are:  (a)  the  signal  strengths  and  S/N 
ratios  for  each  transition  are  high;  (b)  the  signals  strengths  are 
independent  of  source  irradiance  at  either  excitation  line.  The  dis- 
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advantages  of  this  approach  are:  (a)  scatter  of  excitation  radiation  is 
potentially  severe  but  by  measuring  non-resonance  fluorescence  this 
difficulty  can  be  minimized:  (b)  the  post  filter  effect  can  be  severe 
but  it  can  be  minimized  by  measuring  same  fluorescence  wavelength; 
however,  in  the  latter  case,  scatter  can  become  a  problem;  (c)  the 
laser  beams  must  have  the  same  cross-sectional  area;  (d)  the  spectro¬ 
meter  system  must  be  calibrated  at  3-*- i  and  3-j  .  Experimental  results 
for  flame  temperatures  will  be  given  via  this  method. 


Method  3.  Saturation  Two  Line  Method  With  Simultaneous  Pumping 
In  this  method,  both  excitation  transitions  are  pumped  simultan¬ 
eously  in  one  case  and  in  another  case  only  one  excitation  transition 
(1-3  or  2->-3)  is  used.  In  this  case,  the  ratios  of  fluorescence  rad¬ 
iances  is  given  by 


r  3-1 

<^3> 

2-3 

rz 

F3-l 

2-3 


g3  *  8; 

^3 


gl  E12/kT  A31  +  k31 


+  —  e 
®3 


L12 


§1  *  §2  +  §3 
S3 


(5) 


g2  e 

where  all  terms  have  been  defined  above.  Since  k. n  =  k0,  — 

12  21  g1 


-E12/kT 


and  substituting  for  the  g's  (g^  *  g3  =  1;  g2  =  2)  and  solving  for  T 
gives  for  the  case  of  simultaneous  pumping  in  one  case  and  excitation 
at  2-3  in  the  other  case 


T  = 


E^/k 


(6) 


^31  +  k3i 

1  +  4  - — ) 

*21 


Finally,  if  k?1  >>  (A-jj  +  k,^  then, 


if  n2  ^  30  n1 


(7) 


The  same  advantages  as  in  Method  2  applies  here  except  Method  3  is 
also  fairly  insensitive  to  both  scatter  (weak  for  non-resonance  line) 
and  post  filter  effects  (fluorescence  is  at  same  wavelengths)  and  there 
is  no  need  for  spectrometer  calibration.  In  addition,  by  simultaneous 
excitation,  n^,  can  be  evaluated  as  well  as  the  ratio  (A31  +  k^1)/k2j- 
The  major  disadvantage  of  this  method  is  the  need  for  2  laser  beams 
optically  aligned  and  matched  with  respect  to  cross -sectional  area 
(just  as  for  Method  2) ,  and  of  course  the  assumptions  made  above  must 
be  valid.  Experimental  flame  temperature  measurements  by  this  method 
will  be  given  in  this  report. 


Method  4.  Saturation  Two  Line  Method  With  Peak  Detection 

2  ^ 

In  this  method  developed  by  Omenetto  and  Winefordner,  »  it  is 
necessary  to  excite  fluorescence  3-*-l  at  l-*-3  and  a  short  time  later  or 
earlier  3-*-l  at  2-*3.  In  this  case,  the  atomic  system  acts  as  a  2-level 
atom  since  excitation  and  measurement  of  fluorescence  is  done  at  the 
peak  of  the  excitation  profile  prior  to  relaxation  of  the  system  to  a 
3-level  steady  process.  The  temperature,  T,  in  this  case  is  given  by 
the  simplistic  expression 


>  -  r  A 
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The  advantages  of  this  method  are:  (a)  the  k's  and  A's  are  not  needed; 
(b)  the  temperature  obtained  is  independent  of  source  intensity  since 
saturation  must  occur.  The  disadvantages  of  this  method  are:  (a)  the 
spectrometer  system  must  be  calibrated;  (b)  the  laser  beams  must  be 
optically-aligned  and  spatially  matched;  (c)  fast  electronic  detection 
and  impulse  rise  time  of  the  laser  pulses  must  be  achieved  to  allow 
application  of  2-level  steady  state  saturation  theory;  and  (d)  of  course 
saturation  has  to  be  achieved.  This  method  is  the  best  since  it  is  in¬ 
dependent  of  the  approximations  needed  in  the  previous  methods,  except 
that  saturation  is  necessary  and  the  "2-level  peak"  must  be  measured. 

Method  5.  Laser  Excitation  of  Probe-Thermal  Excitation-Emission 
Method  (also  see  Appendix  C) . 

In  this  novel  method,  a  pulsed  laser  is  used  to  excite  an  upper 
level  (say  3)  of  a  probe,  such  as  Tl,  and  then  the  ratio  of  a  thermally 
excited  line  to  a  fluorescence  line  o_r  of  two  thermally  excited  lines 
are  measured.  In  the  first  case,  Method  5a,  the  ratio  of  the  fluores¬ 
cence  at  3-"l  following  excitation  at  2-’-3  to  the  thermally  excited 
emission  at  i-*l  (  where  i  =  4,5,...)  is  taken,  and  so 


1  0 
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In  the  second  case,  Method  5b,  the  ratio  of  2  thermally  excited  lines 
j  +  1  and  i+1  (j  and  i>3)  ,  and  so 
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(Ei3-E.,)/k 

t„./M  ii-iW 

Pf 
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2+3 

l  i"1 

1  2+3 
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In  these  cases,  the  advantages  are:  (a)  only  one  laser  wavelength  is 
needed;  (b )  saturation  is  not  necessary:  (c)  there  is  obviously  no 

need  to  match  laser  beams  in  orientation  and  in  cross  section  area; 

( d)  the  post  filter  effect  is  minimized;  no  scatter  interference;  and 

( e)  since  only  one  laser  excitation  wavelength  is  needed,  temporal  tem¬ 
peratures,  i.e.,  with  only  one  pulse  (^2  ns  -  1  ys,  depending  on  the 
types  of  pulsed  dye  laser)  are  easy  to  obtain  assuming  the  S/N  ratios 
for  the  fluorescence  and  emission  lines  are  sufficient.  The  major  dis¬ 
advantage  is  the  need  to  calibrate  the  spectrometric  system.  The  S/N 
ratio  of  the  emission  lines  will  be  of  the  same  order  as  the  fluores¬ 
cence,  3+1,  due  to  2+3  excitation.  Results  for  flame  temperatures  by 
this  method  will  also  be  given.  The  general  experimental  system  used  for 
flame  temperature  measurements  is  given  in  Figure  2,  and  is  basically 
the  same  as  described  in  the  original  research  proposal. 
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SECTION  III 


THEORY  OF  ABSOLUTE  NUMBER  DENSITY  AND 
QUANTUM  EFFICIENCY  DETERMINATION 

Just  as  for  the  absolute  temperature  measurements,  the  measurement 

.  3 

of  concentrations  (in  species  m  )  requires  the  application  of  2-level 

(peak)  or  3-level  steady  state  theory  (see  Omenetto  and  Wineford- 

~>  3 

ner'1'  and  just  as  for  methods,  2,3,  and  4  for  temperature  measurement 
requires  the  use  of  high  intensity  pulsed  dye  lasers  for  saturation. 

For  a  2-level  atom,  the  fluorescence  radiance  Bp  is  related  [1] 
to  the  continuum  source  spectral  irradiance,  E  ,  by 


where  all  terms  have  been  previously  defined  except  for:  h  =  Planck's 
constant;  £  =  fluorescence  path  length,  and  E^*  is  the  modified  sat- 
uration  spectral  irradiance  (in  J  s  m  4  Hz  )  which  is  given  by 

cA-,.  8-iThv^-,, 

A*  ■  ir2? - z — 21  02) 

v  z  1  *  2 1  c  Y21 

where  B2j  is  the  Einstein  coefficient  of  induced  emission  (transitions 
s  1  (J  m  ^  Hz  *)  *)  and  Y?1  is  the  quantum  efficiency  for  the  fluorescence 
transition  2-*l.  The  saturation  spectral  irradiance,  E^  (i.e.,  the 
source  spectral  irradiance  where  the  fluorescence  radiance  -  50%  of 
the  maximum  possible  fluorescence  under  saturation  conditions) 

is  related  to  E  * 

v 
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(13) 


Ev  '  'in 


By  plotting  Bp  vs  E^  (measuring  absolute  values  of  Bp  and  E^)  ,  E^ 
can  be  determined  and  thus  Y7  ^  can  be  obtained.  Since  Y  ^  p  is  related 
to  radiationless  rate  constants  as  k?p,  such  parameters  can  be 
determined  if  A?p  is  known,  i.e.. 


Y 


21 


21 

~+~~K 
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21 


(14) 


Also  since  k’s  are  related  to  quenching  cross-sections,  one  can  obtain 

quenching  cross-sections  as  long  as  the  composition  and  temperature  of 

the  flame  is  known.  If  either  plots  Bp  vs  E^  or  1/Bp  vs  1/E  ,  it  is 

possible  to  extrapolate  to  E^  -+  00  (or  1/E  +  0)  and  obtain  the  value 

of  B  which  is  given  by 

max 

Bc  =  (4-)  hv17  A71  nT  ( — - )  (15) 

F  v  4  7T  12  21  T  +  pc  J 

max  fcl  s2 


from  which  n^.  can  be  determined  as  long  as  the  g’s,  Vp7,  and  A^p  are 

known  and  SL  and  Br  can  be  absolutely  measured.  Such  measurements 

Fmax 

have  been  used  to  estimate  n^  and  Y  values  of  several  atomic  probes, 
including  Sr,  Ca ,  Na,  In,  etc.  However,  the  values  obtained  have  sys¬ 
tematic  errors  of  as  much  as  a  factor  of  2  to  3  when  compared  with  in¬ 
dependent  measurements  of  and  Y  for  the  same  flames.  In  recent 
weeks,  we’ve  shown  that  the  plots  of  Bp  vs  E^  (or  1/Bp  vs  1/E  )  give 
erroneous  results  for  the  near  plateau  and  early  plateau  of  Bp  and 
E  .  These  conclusions  became  evident  by  systematic  errors  obtained  for 
temperatures  obtained  via  the  use  of  saturation  values  Bp  (Methods  2 
and  3)  obtained  by  plotting  Bp  vs  and  extrapolating  to  the  saturation 
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plateau.  The  reason(s)  for  the  systematic  error  is  (are)  not  known 
to  us  but  are  certainly  related  to  the  complex  effects  of  spectral/ 
spatial/temporal  aspects  of  pulsed  dye  lasers. 


No  attempt  will  be  made  here  to  give  the  corresponding  expressions 

2 

for  3  (or  more)  level  atoms  (one  can  refer  to  Omenetto  and  Winefordner 
for  a  discussion  of  such  cases).  However,  it  is  apparent  that  in 
these  cases,  the  expressions  relating  Bp  to  n^  and  Y ^  also  involve  k's 
and  thus  either  the  k's  must  be  known  or  else  limiting  cases  (such  as 
when  levels  1  and  2  are  very  close  together)  must  be  applied  where  the 
effect  of  the  k'-  is  minimal.  In  addition,  the  same  systematic  errors 
concerning  the  plateau  values  occur  here  also. 


Fortunately,  the  simplest  means  of  obtaining  spatial/temporal  n^, 
values  for  3-level  atoms  is  to  do  it  via  "single"  pulse  saturation. 

For  examples,  if  the  same  conditions  applying  as  for  Method  4  for 
flame  temperatures,  then  n^  can  be  evaluated  from  an  absolute  Bp 
measurement  and  using  the  2-level  Bp  expression  on  the  preceding  page. 
An  innovative  means  of  measuring  n^  for  a  3-level  atom  under  steady 
state  conditions  is  to  saturate  level  3  via  1+3  and  2+3  simultaneously, 
and  using  the  expression  below  (see  Appendices  A  and  B) . 
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<^3> 

2+3 
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One  can  imagine  this  case  as  effectively  as  the  2-level  case  where  levels 
1  and  2  coalesce  into  one  level  of  statistical  weight,  gp  +  g2*  Thus 
here,  the  same  reasoning  as  for  the  2-level  case  applied  and  the  same 
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information  must  he  available  to  determine  n^,.  Unfortunately,  the 
measurement  of  Y,^  is  not  possible  via  the  approach,  and,  as  stated 
above,  it  is  still  necessary  to  know  the  k's  or  to  be  able  to  use 
certain  limiting  cases  where  the  effect  of  the  k's  are  minimized. 
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SECTION  IV 


EXPERIMENTAL  SYSTEM  AND  TECHNIQUES 

All  flames  studies  were  supported  on  a  Meker-type  burner  with  an 
inner  seeded  flame  zone  surrounded  by  an  outer  flame  (unseeded)  of  the 
same  gas  composition.  The  flame  was  surrounded  by  an  inert  gas  sneath 
of  the  same  composition  as  the  major  dilute  gas  species  (i.e.,  Ar  or 
N't).  The  inner  flame  consisted  of  a  rectangular  array  of  32  (4  x  8) 

0  033  in.  dia.  holes  spaced  0.045  in  The  outer  flame  consisted  of 
2  rows  of  holes  surrounding  the  inner  flame.  The  total  flame  zone  had 
8  x  12  holes  spaced  0.045  in  fp.  All  flames  were  run  at  an  unburnt  total 
gas  rise  velocity  of  a-0.4  m  s~*.  The  burner  head  was  water-cooled  with 
25°C  water  and  was  nominally  kept  near  ambient  temperature  during  oper¬ 
ation.  Seed  material  was  introduced  into  the  inner  flame  in  the  form 
of  an  aerosol  produced  from  the  aqueous  solution  of  the  nitrate  or 
chloride  salt  of  Tl,  In  or  Na .  The  aerosol  was  produced  by  either  an 
ultrasonic  nebulizer  (operated  at  ^20  W  power  at  1.4  M  Hz)  or  o>'  a 
Meinhart  type  glass  concentric  nebulizer.  Both  nebulizers  functioned 
adequately.  The  ultrasonic  nebulizer  was  somewhat  less  stable  but 
provided  ^10  X  the  fog  density  to  the  flame  and  was  completely  free 
from  clogging  affects  at  high  salt  concentrations.  Salt  concentrations 
used  ranged  from  0.1  PPM  (s/10^  species  cm  ^  in  the  flame)  to  as  high 
as  1  0,000  PPM  (n,l  0 J  ^  spec  ies  cm  J  in  the  flame). 

The  instrumental  arrangement  is  shown  schematically  in  Figure  3. 

The  N,- laser  used  was  a  Molectron  Corp.  model  UV-14,  in  which  the 
external  power  supply  was  replaced  by  us  with  an  internal  design  power 
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supply  in  order  to  reduce  signal  interfering  radio  frequency  noise 
generated  by  the  ^-laser.  This  ^-pump  was  used  in  pumping  a  Molectron 
Corp.  DL-400  dye  laser  and  a  Lambda  Physik  Model  FL  2000  dye  laser. 

All  single  laser  experiments  employed  the  Molectron  dye  head.  The 
amplifier  stage  of  the  Lambda  Physik  dye  head  was  removed  due  to  the 
insufficient  pumping  power  from  the  .^-laser  after  beam  splitting  the 
UV-14  output.  The  beam  expanding  telescopes  of  both  dye  lasers  were 
detuned  (defocused)  in  order  to  insure  that  a  pseudo-continuum  spectral 
source  profile  would  be  obtained  from  the  laser  output.  The  laser  line 
widths  were  measured  by  scanning  the  laser  wavelength  through  the 
resonance  transition  of  interest  at  low  (linear)  optical  pump  powers 
as  described  by  Omenetto  et  al,  .  Typical  low  power  scans  for 
Agx  ^->3  and  Xex  2+3  an  are  shown  in  Figure  4.  The  beam  was  then 
passed  through  a  series  of  circular  spatial  filters  in  order  to  achieve 
a  spatially  uniform  intensity  distribution. 

The  detection  system  consist  of  a  wide  bandwidth  monochromator 
(Jobin  Yvon  Model  H-10)  employed  as  the  wavelength  selection  device. 

The  signal  was  collected  by  a  Hamamatsu  R928  photomultiplier  take  with 
a  specially-wired  base  for  high  speed  transient  work.  Output  signals 
were  found  to  be  linear  with  light  intensity  up  to  40  mA  peak  current 
with  no  temporal  distortion;  beyond  this  current,  the  peak  would  become 
temporally  distorted,  and  only  the  time  integrated  signal  was  found  to 
be  linear  with  light  intensity.  Improvements  on  this  response  could 
extend  the  linear  peak  current  region  to  ^1  A,  and  thereby  extend  the 
expected  S/N  ratio  on  single  pulse  measurements.  The  PMT  was  operated 
with  an  uneven  dynode  chain  at  -1400  V  potential  (cathode  to  anode). 
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,  .  3-2/1-3 

both  cases,  (see  Appendix  D) . 


The  measured  rise  time  of  a  250  ps  input  pulse  was  found  to  be 
x800  ps  (for  95%  signal).  This  signal  was  fed  into  a  Tektronix  Corp. 
Model  151  sampling  oscilloscope  with  350  ps  rise  time.  Our  system's 
overall  temporal  response  (for  951  signal  level)  is  M . 2  ns  which  was 
sufficient  for  the  temporal  studies  and  data  evaluation  contained 
in  this  report.  In  order  to  minimize  jitter  problems  with  the  Mole- 
ctron  N-,  laser,  an  optical  triggering  arrangement  was  employed.  This 
triggering  scheme  eliminated  the  possibility  of  using  the  PAR  Boxcar 
detection  systems  in  our  laboratory  due  to  their  requirement  of 
75  ns  pretrigger  pulse  (all  attempts  at  signal  delay  of  75  ns  caused 
distortion  of  the  temporal  waveform). 


All  electronic  detection  systems  were  connected  to  a  salt-bed 
earth-ground  in  order  to  minimize  R.F.I.  problems.  The  noise  limiting 
our  signal  detectability  in  this  system  was  found  to  be  noise  in  the 
electronic  measurement  system  for  the  non-resonance  cases  and  scatter 
noise  in  the  resonance  cases.  The  detectability  limit  for  T1  with 

*7 

system  was  found  to  be  ^0.5  part-per-trillion  in  solution  or  %5.  x  10 

-  3  3 

atoms  cm  in  the  flame  using  a  detection  volume  of  0.5  x  1.5  x  3  mm  ; 

this  volume  corresponds  ^10-100  atoms  as  the  absolute  detectability 

limited  for  the  system  described.  All  measurements  were  made  with  an 

averaging  time  constant  of  0.5  s  and  a  repetition  rate  of  20  Hz. 

For  the  linear  and  thermally  assisted  temperature  measurements,  the 

species  density  was  %5  x  10  cm"  ,  while  ^5  x  10  cm  number  density 


was  employed  for  the  saturation  measurements.  (The  1.68  eV  thermally- 
assisted  measurement  re< 
reasonable  S/N  ratios) . 
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assisted  measurement  required  ^5  x  10  cm  densities  to  obtain 
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No  "single  pulse"  temperatures  are  given  (we  did  not  have  avail¬ 
able  the  2  dye  lasers  after  we  had  evaluated  the  meaning  of  the  tem¬ 
poral  photomultiplier  scans  for  the  saturation  temperature  methods). 


However,  we  did  measure  several  "single  pulse"  fluorescence  signals 
by  METHOD  # 2  using  In  as  a  probe  (this  probe  was  later  shown  to  give 
too  large  random  temperature  errors  due  to  the  small  separation  of 
levels  1  and  2)  and  obtained  reasonably  good  fluorescence  signals 
(the  uncertainty  in  the  fluorescence  signals  would  correspond  to  a  tem¬ 
perature  error  of  <100  K  with  T1  as  a  probe).  We  should  also  point 
out  that  all  flame  temperatures  correspond  to  only  about  10  laser 
pulses  for  each  fluorescence  signal  indicating  the  rather  excellent 
signal-to-noise  ratios  with  the  saturation  methods. 


21 


The  experimentally  measured  flame  temperatures  obtained  by 
Methods  1-5  and  a  comparison  with  previously  measured  values  by  the 
Utrecht  group  and  with  the  sodium  line  reversal  method  in  one  labora¬ 
tory  are  given  in  TABLE  I.  The  footnotes  in  Table  I  indicate  the  source 
of  random  errors.  In  TABLES  II,  III,  and  IV,  the  influence  of  certain 
experimental  parameters  on  the  flame  temperatures  measured  by  Methods 
#1  and  #5  are  given.  The  correlation  of  temperature  values  measured 
by  Methods  #1-5  with  previously  measured  values  is  generally  very  good. 
There  are  obvious  sources  of  systematic  errors  even  between  the  Utrecht 
values  as  well  as  between  the  values  obtained  by  us  with  the  5  methods. 

The  estimated  random  temperature  errors  in  all  of  the  present 
studies  were  made  conservatively,  i.e.,  using  3/5  peak-to-peak  noise 
rather  than  1/5  peak-to-peak  noise  as  is  recommended  by  most.  It 
should  be  possible  to  reduce  the  temperature  errors  from  ±  30  K  to 
+  10  K  (or  less)  in  the  saturation  method  by  using  a  wider  gate  (2  ns 
vs  350  ps)  and  extracting  higher  peak  currents  from  the  photomultiplier 
detector.-  In  the  thermally  assisted  temperatures  (METHOD  #5)  and  the 
linear  method  (METHOD  1),  the  random  temperature  errors  due  to  S/N 
and  calibration  can  also  be  reduced  to  nearly  +10  K.  It  should  be 
stressed  that  all  flame  temperature  measurements  reported  here  are 
preliminary  ones  and  can  be  further  refined,  particularly  with  the 
present  knowledge  if  the  mechanistic  aspects  resulting  from  the  present 
studies  and  the  information  to  be  obtained  in  future  studies  (refer 
to  CONCLUSIONS) . 


TABLE  II.  EXPERIMENTAL  CONDITIONS  AND  INFLUENCE  OF  SOURCE  IRRADIANCE  ON  LINEAR  TWO-LINE 
(METHOD  #1),  FLAME  TEMPERATURES.  (Also  refer  to  Figure  3). 
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TABLE  IV.  INFLUENCE  OF  LASER  IRRADI ANCF.  UPON  THERMAL  EXCITATION  EMISSION  (METHOD 
TEMPERATURES.  (Also  refer  to  Figure  4). 


In  Figures  5  and  6,  typical  Bp-E^  plots  are  given  for  several 
fluorescence  transitions  T1  in  several  H^-based  flames.  In  Method 
I,  Bp  values  are  measured  for  any  source  spectral  irradiance  on  the 
linear  portion  of  the  plots.  In  Methods  2  and  3,  Bp  values  must  be 
taken  on  the  3-level  steady  state  plateau  region,  and  in  Method  4 
Bp  values  must  be  taken  on  the  2-level  steady  state  plateau  region. 

In  Method  5,  Bp  values  can  be  taken  in  any  portion  of  the  Bp-Ey  plots 


In  the  linear  two  line  fluorescence  method  (Method  I),  it  is 
critical  to  have  an  accurate  and  precise  calibration  of  source  irradiances 
and  fluorescence  radiances;  the  former  requires  a  careful  calibration 
of  the  laser  spectral  bandpass  for  both  excitation  transitions.  In 
addition,  it  is  necessary  to  assure  oneself  that  the  Bp  values  corres¬ 
pond  to  the  linear  part  of  the  Bp.  vs  E  curves  (see  TABLE  II  and 
Figure  3).  The  major  random  errors  occur  due  to  the  poorer  signal-to- 
noise  ratios  (as  compared  with  Methods  based  on  Saturation)  unless 
higher  number  densities  are  employed  (with  the  concomitant  problems 
associated  with  scatter,  post  filter,  etc.),  calibration  errors,  and 
laser  spectral  bandpass  errors.  One  must  not  estimate  ^^jaser  since 
an  error  of  0.004  nm  can  result  in  a  flame  temperature  error  of 
%200  K,  and  since  this  approach  requires  a  pseudo -cont inuum ,  it  is 
necessary  to  detune  to  spectrally  broaden  the  laser  output.  Because 
of  the  need  to  match  approximately  laser  beam  geometries  and  the  smaller 
S/N-ratios,  this  approach  is  less  suitable  for  measuring  temporal 
(<"  1  us  time  period)  flame  temperatures.  In  Methods  2-4,  it  is  necessary 
to  know  whether  the  measurement  with  a  suitable  detector  gate  corresponds 
to  either  the  "peak"  (METHOD  4)  or  to  the  "steady  state"  (METHODS  2 
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and  3).  In  Figures  7a  and  7b,  the  temporal  scans  of  the  fluores¬ 


cence  response  are  given  for  the  2  flames  (I^/C^/Ar  and  It 

is  apparent  that  in  the  less  quenching  Ar-diluted  flame  that  2  regions 

corresponding  to  the  "peak"  and  "steady  state"  cases  exist,  whereas 

in  the  more  quenching  ^-diluted  flame  only  one  region  occurs  partly 

due  to  the  inadequate  temporal  resolution  of  the  "peak"  and  "steady 

state"  regions  via  our  detection  system.  Flame  temperatures  obtained 

from  the  scans  in  Figures  7a  and  7b  are  given  in  TABLE  V.  In  METHODS 

2,3,  and  4,  beam  matching  is  also  necessary  as  the  necessity  of  a 

pseudo - cont inuum  source.  The  exact  value  of  AA,  is  not  necessary 

r  laser  7 

as  long  as  the  dye  laser  approximates  a  pseudo-continuum  with  respect 
to  the  absorption  line.  With  commercial  dye  lasers,  it  is  necessary 

o 

to  detune  (by  0.1  or  0.2  A)  in  order  to  approximate  a  pseudo-continuum, 
but  in  so  doing,  the  laser  irradiance  is  reduced,  in  some  cases  suffi¬ 
ciently  to  no  longer  saturate  the  excitation  transition.  Therefore, 

the  excimer  laser-dye  laser  system  is  needed  (currently  cur  NS-dye  laser 

-  2 

system  is  capable  of  peak  spectral  irradiances  of  about  1-20  MW  cm  nm 
in  the  fundamental  region  and  about  1/10  to  1/100  of  that  in  the  doubled 
region;  the  excimer  pumped  dye  laser  will  have  spectral  irradiances 
of  the  order  of  100X  greater  than  the  present  N9-dye  laser  system  and 
should  be  adequate  to  saturate  all  visible  and  most  ultraviolet 
resonance  transitions).  In  addition,  in  METHODS  2,3,  and  4,  scatter 
and  post  filter  errors  must  be  of  concern.  The  best  excitation  and 


fluorescence  transitions  involve  B„  and 


since  scatter  and 


post  filter  errors  are  minimized  in  the  former  and  are  small  in  the 
latter  since  the  fluorescence  signals  are  rather  large  and  there  is 


1 


i 


no  need  to  know  gA-values*. 

Method  5  has  the  greatest  number  of  desirable  characteristics  for 
single  pulse  temperature  measurements  since  there  is  no  need  to  beam 
match  or  saturate.  However,  there  is  need  to  calibrate  the  spectrometric 
system.  The  greatest  difficulty  with  this  method  is  the  need  (like 
in  Method  1)  to  know  gA-values*  which  are  oftentimes  not  accurately 
known.  Poorly  defined  gA-values  lead  to  systematic  temperature  errors. 

It  is  possible  in  this  method  that  both  peak  and  integrated  pulse 
measurements  will  result  in  accurate  temperatures.  Studies  on  this 
will  be  done.  Also,  the  most  critical  study  to  be  performed  involves 
the  need  to  verify  thermal  population  of  the  thermally-excited  levels 
(4,5,...)  in  the  probe,  i.e.,  experimental  measurements  of  collisional 
rate  constants  k^  when  j  =  4,5,...  Such  studies  will  be  completed 
in  the  near  future. 

Finally,  in  Figure  8,  indirect  evidence  is  given  indicating  the 
possibility  of  thermal  population  of  level  4.  The  fluorescence  pulse 
precedes  the  thermal  population  of  level  4.  The  fluorescence  pulse 
precedes  the  thermally  excited  emission  pulse  by  only  ^1  ns  which 
lends  support  to  equilibration  of  the  excited  levels  within  a  few 
nanoseconds . 

We  have  also  determined  that  plots  of  Bp  vs  resulting  in 
"saturation"  plateaus  give  erroneous  values  of  flame  temperatures  T's, 

*A11  gA  values  taken  from  C.H.  Corliss  and  W.R.  Bocman,  NBS  Monograph 
53,  U.S.  Dept,  of  Commerce,  1962. 
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and  concentrations,  n^'s;  the  saturation  plateaus  resulting  from 
extrapolation  techniques  (i.e.,  dn^/dt  f  0,  dn2/dt  f  0  during  the  N^-dye 
I  laser  excitation  pulses)  and  from  sufficiently  high  laser  spectral 

irradiances,  have  apparent  systematic  errors.  Therefore,  previous 
methods  based  on  measurement  of  n^'s  of  OH,  CN,  CH,  etc.,  from  satura¬ 
tion  plateaus  easily  by  a  factor  of  2-3X  in  error. 


SECTION  V 


CONCLUSIONS 

/ 


f 


1.  Theoretical  development  of  4  new  methods  of  nteasuring  spatial 
flame  temperatures  and  comparison  of  these  'A  methods  with  the 
linear  two  line  fluorescence  method  as  well  as  the  line  rever¬ 
sal  method.  The  5^  fluorescence  based  methods  are: 

Method  1.  Linear  Two  Line  Method 

Method  2.  Saturation  Two  Line  Method  with  Sequential  Pumping 

Method  5.  Saturation  Two  Line  Method  With  Simultaneous  Pumping 

Method  4.  Saturation  Two  Line  Method  With  Peak  Detection 

Method  5.  Laser  Excitation  of  Probe-Thermal  Excitation  -  Emission 

Method . 

2.  Experimental  verification  of  the  above  theoretical  approaches  to 
measurement  of  flame  temperatures  of  H?-based  flames.  Problems, 
errors,  and  concerns  associated  with  each  method  are  evaluated. 
Extensions  of  the  5  methods  to  more  quenching  flames  is  given; 
basically  METHOD  1  and  5  can  be  extended  to  hydrocarbon  flames 
with  mainly  a  problem  in  the  reduction  of  fluorescence  signal 
levels  when  operating  on  the  linear  portion  of  the  Bp  vs  E^ 
curves.  METHODS  2 ,  .3 ,  and  4  require  saturation  to  be  achieved; 
conventional  N9-dye  laser  systems  are  incapable  of  saturation  of 
excitation  transitions  below  %500  nm.  Excimer-dye  laser  systems 
possess  sufficient  spectral  irradiance  to  saturation  over  UV 
excitation  processes  in  most  quenching  flames. 
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3.  Theoretical  development  of  2  new  approaches  for  measuring  the 
total  number  density,  n^,,  of  species  in  flames.  The  2  methods 
are  based  on  the  "peak"  2-level  method  and  the  simultaneous 
saturation  of  an  upper  level  by  2  laser  excitation  processes. 
Problems  associated  with  the  conventional  1/Bp  vs  1/E^  plot 

to  determine  n^,  and  Y  (quantum  efficiency)  are  indicated.  The 
"saturation"  plateau  fluorescence  radiances  are  found  to  give 
systematic  errors  on  both  the  n^,  and  T  values,  and  so  it  is 
recommended  that  saturation  fluorescence  signals  be  measured 
via  either  the  "peak"  or  "steady  state"  methods  (see  Figures 
7a  and  7b) . 

4.  Other  studies  completed  include: 

a.  Theory  and  Experimental  evaluation  of  Saturation  broadening. 

b.  Experimental  evaluation  of  n^-values. 

5.  Manuscripts  to  be  forwarded  (when  completed) 

a.  Theory  of  Fluorescence  Temperature  Methods 

b.  Experimental  Verification  of  Fluorescence  Temperature  Methods 

c.  Theory  of  Fluorescence  Number  Density  Methods 

d.  Saturation  Broadening 

e.  Experimental  Evaluation  of  n^ 

Summary  of  Studies  To  Be  Completed 

1.  Flame  temperature  of  hydrocarbon  based  flames  (gaseous  and  liquid 
fuel)  will  be  completed  as  soon  as  the  new  excimer  laser  is 
available.  Temporal  (<  1  ps)  based  flame  temperatures  will  also 
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be  done  when  the  excimer  laser  is  available.  The  excimer  laser 
will  allow  satisfactory  pumping  of  2  dye  lasers  for  the  satura¬ 
tion-based  methods.  "Temporal"  temperatures  were  not  obtained 
in  the  present  work  (other  than  proof  that  the  fluorescence 
signals  were  sufficient  to  obtain  accurate  T's)  because: 
fa)  2  dye  lasers  were  available  only  for  3  months  of  the 
contractual  time  and  (b)  only  after  return  of  the  Lambda  Physik 
dye  laser  to  INRAD  did  we  realize  the  basis  of  METHOD  5  and  the 
meaning  of  the  transient  fluorescence  signals  (see  Figures  6 
and  and  discussion  concerning  time). 

2.  A  more  thorough  study  of  the  transient  fluorescence  signal  is 
still  needed  to:  (a)  prove  thermal  population  of  the  thermally 
excited  levels  in  METHOD  5,  and  (b)  determine  whether  the  "peak" 
and  3-level  stead)-  state  portions  of  the  fluorescence  transient 
signals  can  be  resolved  via  a  faster  detection  system  than  we 
now  possess.  For  example,  we  plan  to  probe  level  2  of  T1  by 
pumping  l-*-3  with  a  N-,  laser  and  using  a  Chromatix  CMX-4  to 
monitor  level  2  by  pumping  2-*4  (or  other  level). 

3.  Finally,  we  hope  to  transfer  over  methods  to  native  molecular 
species  in  flame,  e.g.,  OH,  CH,  CN,  C7 ,  etc. 
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SECTION  VI 


MISCELLANEOUS  ITEMS 


PERSONNEL  ON  PROJECT 

The  present  project  was  performed  by  one  post  doctoral  associate 
(Dr.  Nicolo  Omenetto  now  at  CCR,  European  Community  Center,  Stabilimento 
di  Ispra,  Ispra,  ITALY)  for  12  months  paid  directly  from  contract  funds 
and  two  full  time  graduate  assistants  (for  18  months)  paid  primarily 
via  funds  from  other  contracts  (John  Bradshaw  and  James  Bower) .  In 
addition,  Dr.  Benjamin  Smith;  Dr.  Georgio  Zizak  -  Centro  di  Studio  Per 
Ricerche  Sulla  Propulsion  E  Sull'  Energetica,  Milan,  ITALY  ;  and  Mr. 

John  Horvath  have  spent  4  months  and  10  months,  respectively,  on  pro¬ 
jects  involved  in  this  contract. 

NON -EXPENDABLE  EQUIPMENT  PURCHASED 

Equipment  (non-expendable)  purchased  on  this  contract  includes: 

(a)  Lambda  Physik  dye  laser  for  $15,000  from  Inrad,  Inc., 

181  Legrand  Ave.,  Northwale,  NJ  07647;  (b )  Portable  Gas  (Nitrogen) 
Supply  System  for  liquid  nitrogen  for  $1125  from  Union  Carbide  Corp., 
Linde  Division,  4800  West  16th  Street,  Indianapolis,  IN  46244;  (c) 

Kinetic  Systems  Gated  Connector,  Unibus  Cable,  Extender  Card,  Camac 
Crate  Controller  and  Camac  Crate  with  power  supply  from  Kinetic  Systems 
11  Mary  Knoll  Drive,  Lockport,  IL  60441;  (d)  LeCroy's  12  Channel  A/D 

Converter  with  30  ns  gate  width  to  fit  in  the  CAMAC  Crate  from  LeCroy 
Research  Systems,  700  South  Main  St.,  Spring  Valley,  NY  10977;  (e ) 

RT-11  Operating  System  with  Fortran  Extensions,  Scientific  and  Lab 
Subroutine,  Basic  Extensions,  Lab  Applications  and  6  Disk  Cartridges; 
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and  (e)  several  small  items  ($220  or  less)  including  a  micro  pump 
and  a  gas  regulator. 

The  Lambda- Phys ik  dye  laser  had  arrived  about  9  months  after  the 
start  of  the  contract  and  the  dye  laser  scan  unit  more  than  1  year 
after  the  start  of  the  contract.  The  combined  system  was  available 
for  the  dual  dye  laser  excitation  measurements  (Molectron  DL-400  was 
other  dye  laser)  for  less  than  3  months  out  of  the  total  contract 
period.  To  make  matters  worse,  the  dye  laser  system  never  met  specifi¬ 
cations  concerning  spectral  resolution  step  size,  wavelength  accuracy 
and  repeatability,  wavelength  range,  and  several  other  less  serious 
aspects.  Finally,  in  November  1979,  the  Lambda - Phys ik  dye  laser  was 
returned  to  the  vendor  (Inrad)  and  the  entire  dye  laser  system  was 
rebid  (within  the  past  2  weeks  the  new  dye  laser  -  Molectron  DL-14) 
has  arrived. 

All  other  major  items  were  required  for  temporal  fluorescence 
measurements  and  computer  interfacing. 


REFERENCES 


f 

A 


* 


1 

i 


1.  N.  Omcnctto,  F.  Bcnetti,  and  G.  Rossi,  Spectrochim.  Acta,  27-b 
453  (1972). 

2.  N.  Omcnctto,  R.F.  Browner,  J.D.  Winefordner,  G.  Rossi,  and  P. 
Bcnctti,  Ana  1 .  Chem. ,  44 ,  1683  (1972). 

3.  P.  Bcnetti,  N.  Omcnetto,  and  G.  Rossi,  La  Termtecnica,  29,  86 
(1975). 

4.  H.  Haraguchi,  B.  Smith,  S.  Weeks,  D.J.  Johnson,  and  J.D. 
Winefordner,  Appl icd  Spectroscopy ,  3 1 ,  156  (1977). 

5.  H.  Haraguchi  and  J.D.  Winefordner,  Applied  Spectroscopy,  31,  195 
(1977)/ 

6.  H.  Haraguchi  and  J.D.  Winefordner,  Applied  Spectroscopy,  31,  330 
(1977)  . 

7.  J.  Bradshaw,  J.  Bower,  S.  Weeks,  K.  Fujiwara,  N.  Omenetto,  H. 
Haraguchi,  and  J.D.  Winefordner,  10th  Materials  Research  Symposium 
on  Characterization  of  High  Temperature  Vapors  and  Gases,  NBS^ 

Ci'a  ithershurg  ,  Maryland  ,  1978  . 

8.  M.  Lapp  and  C.M.  Penney,  Editors,  "Laser  Raman  Gas  Diagnostics", 

P 1 enum  Press,  1974. 

9.  R.M.  Measures,  J.  Appl .  Phys. ,  39,  5232  (1968). 

10.  N.  Omenetto  and  J.D.  Winefordner,  Chapter  4  in  "Analytical  Laser 
Spectroscopy",  N.  Omenetto,  Editor,  Wiley.  1979. 

11.  N.  Omcnetto  and  J.D.  Winefordner,  Progress  in  Analytical  Atomic 
Spectroscopy,  Vol .  2  (1,2),  Pergamon,  1979. 

12.  G.D.  Boutilier,  M.B.  Blackburn,  J.M.  Mermet,  S.J.  Weeks,  H. 
Haraguchi,  J.D.  Winefordner,  and  N .  Omenetto,  Applied  Optics, 

1_7  ,  2  291  (1  978)  . 

13.  N .  Omcnctto,  J.N  Bower,  J.D.  Bradshaw,  C.A.  Van  D  i  i  k ,  and  J.D. 
Winefordner,  J.  Quant .  Spec t  rose .  R  a  d  i  a  t .  Trans  fer ,  in  press. 


i 


4  4 


APPEND  I  X  A 


(APPL.  OPTICS,  17,  2291  (1978)). 


STEADY  STATE  ATOMTC  FLUORESCENCE  RADIANCE 


EXPRESSIONS  FOR  CONTINUUM  EXCITATION* 


C.  D.  Routilier 
M.  B.  Blackburn 
J.  M.  MernefC 
S.  j.  Meeks 

£ 

H,  HsraPuchi 

J.  D.  Winefordner^ 

Department  of  Chemistry 
University  cf  Florida 
Gainesville,  FE  32611 


N,  Omenetto 

Institute  of  Inorganic 
and  General  Chemistry 
University  of  Pavia 
Vi  ale  Taramelli  12 
27100  Pavia 
Italy 


Work  supported  bj'  Af-AFOSR-F4',i6  20-76-C-0005. 

G.  1).  B.  wislie.;  to  acknowledge  the  support  of  a  fellowship  sponsored 
by  the  Procter  (■  Gamble  Co. 

c.  Present  address:  Lab.  de  Phy riccchimie  Ind. ,  Bat  401,  INSA,  69621 

Vil  lourfcamio  CF.DEX,  France 

d.  Present  address:  Analytical  Chemistry  Division,  National  Bureau  of 

Standards,  Washington,  D.C,  20234 

e.  Present  address:  Department  of  Chemistry  and  Physics,  National  Institute 

for  Environment,  R.O.  Yarabc,  Ibarakl  300-21,  Japan 


f.  Avithor  to  whom  rrprir. t  is  should  he  sent. 


5 


ABSTRACT 

Several  general  expressions  for  the  atomic  fluorescence  radiance 
excited  by  a  laser  source  are  presented  for  both  2-level  and  3-level  systems. 
The  laser  has  been  approximated  to  a  spectral  continuum  since  its  bandwidth 
is  assumed  to  be  much  larger  than  the  absorption  profile.  Two  types  of 
3-level  systems  are  considered:  i)  systems  where  2  excited  states  are 
close  together,  both  coupled  radiatively  with  the  ground  state  (Na,  K,  Cs, 
etc.);  and  ii)  systems  where  one  level  is  close  to  the  ground  state,  both 
being  radiatively  connected  with  a  third  excited  state  (Ti,  In,  Ga,  etc). 
Expressions  are  given  for  both  linear  and  saturating  conditions  and  the  use 
of  such  expressions  for  absolute  species  measurement  is  discussed.  In  our 
treatment,  we  have  made  the  assumption  that  the  rate  equations  approach  is 
valid. 
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INTRODUCTION 


The  at.  tractive  feasibility  of  laser  excited  fluorescence  spectroscopy 

for  plasma  diagnostics  has  been  described  early  in  1968  by  Measures'*'  and 

recently  exploited  by  a  number  of  researchers  for  measuring  spatial 

species  concentration,  temperature,  and  other  physical  parameters  of 
2_6 

interest.  The  laser  systems  which  can  he  utilized  in  these  studies  are 

the  N^-pumped  dye  laser,  f lashlar.p  pumped  dye  laser,  and  the  cw  Ar  ion- 

purnped  dye  laser.  Very  often  the  laser  bandwidth  is  much  larger  than  the 

absorption  profile,  irrespective  of  the  type  of  broadening  occurring  in  the 

plasma  investigated  and  therefore  it  can  be  approximated  to  a  spectral 

A 

continuum.  For  this  excitation,  Omenetto,  et  al. ,  have  evaluated  the 
2-level  atom  fluorescence  radiance  expressions.  Omenetto  and  Winef  ordner'* 
have  discussed  both  the  2-level  and  3-level  cases  and  Daily^’^  and 

g 

D'Olivares  have  also  evaluated  both  the  2-  and  3- level  cases.  Since  the 
expressions  given,  apart  from  the  different  terminology,  are  of  limited 
applicability,  the  aim  of  this  note  is  to  extend  the  derivation  to  other 
typical  3-level  atomic  systems  and  to  indicate  the  use  of  such  expressions 
for  absolute  species  measurement.  The  effect  of  saturation  with  high 
intensity  laser  sources  can  then  be  consistently  compared  assuming  the  atom 
approximates  either  a  2-level  or  a  3-level  system. 

It  is  worth  stressing  that  the  rate  equations  approach  is  considered 
to  be  valid  here.  Finally,  the  expressions  given  are  only  applicable  in 
the  case  of  low  optical  densities  (k^i  <  0.05). 

atomic  fluorescence  radiance  expressions 

The  atomic  fluorescence  radiance  for  any  atom  (2  or  multi-level)  is 
given^  by 
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b  =  [—)  y  E  B.  ni  [i  -  —  -E1) 

1H  Pjk  viml  c  )  lm  1  K  n #. 


where: 


B_  **  atomic  fluorescence  radiance  for  j-*k  radiative  deactivation 

Vic  _x 

1-mu  process  and  for  i-*-m  radiative  excitation  process,  J  s 

-2  -1 
m  sr  ; 

l  *=  fluorescence  path  length  in  direction  of  observation,  m; 

Y  =  fluorescence  power  efficiency,  J  s  ^  fluoresced/J  s  ^  absorbed 
Pjk 

(where  j  is  upper  level  of  fluorescence  atid  k  is  lower  level); 

Y  =  f lvorcscence  quantum  efficiency  (see  Appendix  A),  dimension- 
3  k 

less  (Y  =  Y  hv  /hv  ) ; 

Pjk  jk  jk  im 

E  =  spectral  irradiance  of  continuum  source  at  frequency  v.  of 
Vira  m 

excitation  (£  is  lower  level  from  which  transition  begins  and 

-1  -2  -1 

m  is  upper  level  at  which  the  transition  ends) ,  J  s  m  Hz 

hv  =  energy  of  exciting  photon,  J; 
c  =  speed  of  light,  tr.s-*; 

B  =  Einstein  coefficient  of  induced  absorption,  i.e.,  absorption 

transitions  per  spectral  energy  density  per  absorbing  species, 

3  -1  -1  u 

m  J  s  Hz; 

n^  =  density  (concentration)  of  species  (absorbers)  in  lower 
level  £,  ro  ; 


3i  or  m 


=  statistical  weights  of  £  or  in  levels,  dimensionless ;  and 


n  =  density  (concentration)  of  upper  level  of  excitation,  m 
m 

An  alternate  expression  for  B  is 


v  -M  »■ 

1  .  ,  law  j 
jrl:  V 


hv  n. 
Jk  jk  j 
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For  a  2  level  atom,  where  the  levels  are  designated  ground,  1,  and, 


1st  excited,  2,  the  terms  in  equation  ]  are:  Y  ,  E  ,  n 
hv  ,  B12’  aud  S° 


n  ,  n  *=  n_, 
1  m  2 


21 


£]  hi  h, 


12' 


hv 


12 


gl  n2 

B.,  [1  -  ~+] 

12  1  82  nl 

3,4 


(2) 


Assuming  steady  state  excitation/de-excitation,  assuming  negligible  popula¬ 
tion  in  other  levels  (n^,  =  +  n^)  *  and  assuming  negligible  thermal  excitation 

E 


g2 

"2  '  &1  nl 


12 


E  +  E*‘ 

V  V 

12  12 


nl  = 


Ht 


1  +  — 


12 


+  E 


12 


* — 
v12-) 


(3) 


where: 


c  A 


21  .  -1  -2  „  -1  . 
,  J  s  m  Hz  ;  and 


V12  Y21B21 

=  Einstein  coefficient  of  induced  emission  (B2lg2  *  B12gl^ * 

3  t-1  -1  „ 

m  J  s  Hz. 

* 

The  term  E  is  a  modified  saturation  spectral  irradiance  and  is  related  to 

vi  2 

*•  ^  8  9 

the  one  given  by  Omenetto  et  al.,‘  by  D'Olivares,  and  by  Daily  by 


=  E" 


f®2  S-^1 


12 


12 


(4) 


The  previously  used  saturation  spectral  irradiai.ee,  F.  •,  will  not  be  used 

12 

here  despite  its  obvious  physical  usefulness  (i.e.,  E®  is  that  source 

12 

spectral  irradiance  where  the  fluorescence  signal  is  50%  of  the  maximum 

possible  value)  because  E‘  is  simply  related  to  E  (evaluated  from 

V12  V12 
fundamental  constants,  c,  and 


measureable  parameters,  A,^,  B.^*  ^2j^ 
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w  v. • 


only  for  the  2-level  system. 


The  general  B  expression  in  terms  of  n_,  for  a  2-lcvel  system  is 
*21  T 


21 


*  \ 

hv  .) 

4?  v21  Ev19 

12 

c 

)  12 

t  4 

glnT 

*12  [  "(g.+g.HE  7e*~  ')"  +  iT 

A  z  12  12  l 


)  (5) 


For  a  low  intensity  source  (E  <<  E  ) 

V12  V12 


B_  “  T~  A,,  hv.„ 

F21(1o)  471  21  12 


Ei 


12 


12 


2  ,  l  v  _ 

g1  °r  4  it  Y21B12 


hv 


12 


Ev  V 
12  1 


For  a  high  intensity  source  (E  >>  E  ) 

12  V12 


B„  "  7T  A91  hvio 

F21(hi)  4  21  12 


f_ ?U 

1.8^  "h  g2J 


(7) 


6  7  2-5 

Both  expressions  are  consistent  vjith  those  of  Daily  *  and  Omenetto. 


Three  Level  System — General  Expressions 

There  are  two  major  types  of  3- level  atomic  systems:  (I)  systems 
where  the  levels  2  and  3  are  close  together  and  radiative  transitions 
between  levels  2  and  3  are  not  allowed,  e.g.,  Na,  K,  Cs ,  etc.;  and  (II) 
systems  where  levels  1  and  2  are  close  together  and  radiative  transitions 
between  1  and  2  are  not  allowed,  e.g.,  Tl,  In,  Ga,  etc.  In  Appendix  A,  the 
general  steady  state  approach  for  estimating  the  ratios,  n^/n^,  and 

n2^nl  £*Aev1,  Below  .only  specific  expressions  for  these  ratios  will  be 
given.  The  specific  cxcitation-deexcitution  routes  are  given  in  parentheses. 

Case  TA(l-»3).  In  this  case,  the  non- rad i ative  excitation  rate  constants 
k^2  and  are  assumed  negligible  and  the  radiative  rate  constant,  =  0, 

Therefore,  thermal  excitation  is  negligible,  and  assuming  radiative  excitation 
of  level  3  from  level  1,  the  ratios  of  n^/n^  and  n2^nl  art? 


(6) 
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V 


n3  _ 

f3 

f  E 

V13 

nl  " 

81 

E  +  K* 

V13  V13' 

n2  _ 

£3 

r  E 

V13 

k 

32 

nl  ’ 

81 

E  +  E*  ~ 

V13  V13' 

A21+k21+k23. 

(8) 


(9) 


where: 


13 


'13 


excitation  spectral  irradiance  at  exciting  frequency  v^. 


.  -1  -2 
J  s  m  ; 


spectral  irradiance  needed  to  make  the  density  of  n3  50%  of 

.  .  -1  -2 
n1  assuming  **  83*  J  s  m  ; 


A^  =  Einstein  coefficient  of  radiation  emission  for  radiative 


deactivation  process,  2-*l,  s  ^ ; 


-1 


k's  =  non-radiative  rate  constants  for  processes  shown,  s  ;  and 

g's  =  statistical  weights  of  levels  shown,  dimensionless. 

* 

The  term  E  is  defined  by 
V13 


cA 


31 


CS3A31  , 


V13  Y31E31  s1Y3IB13 


(10) 


where: 


B^  =  Einstein  coefficient  of  induced  emission  for  process  3+1, 

3  t-1  ~l  t,  . 

tn  J  s  Hz; 

=  Einstein  coefficient  of  induced  (radiative)  excitation  for 


u,  3  -1  -1  u 
process,  l->-3,  m  J  s  Hz; 


10 


and  ^31’  the  fluorescence  quantum  efficiency  for  the  radiative  process  3-+-1 


is  given  by 
Y31 


*31 


k3?k23 

A  +  Yr  +  1-  M  - — — — - - - 

31  j]  ‘32  (A3i-H:3i+U32)(A?Jvk?i+kr3) 


(ID 
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n..  I'.^ni 


whore  the  term  in  brackets  in  the  denominator  accounts  for  the  increase  in 
quantum  efficiency  via  the  cyclic  factor^of  crossover  from  3-0  and  re¬ 
crossover  from  2  0. 

Case  I  Aa(l-Q,  30)  .  Solving  for  n^  in  terms  of  n^,  (n^,  =  n]^  +  n2  +  n3^ 
utilizing  equations  8  and  9  and  substituting  for  n^  into  equation  la  gives 

1 


'  3-*-l 
1-0 


l 


A31  hv13  n,r 


1  + 


32 


Bi 


A21+k21+k23 


1  + 


13 


iy 


(12) 


Case  I  Ah  (1-0,  3-0,  20).  Solving  for  n^  in  terms  of  n^,  utilizing  equations 

10 


8  and  9  and  substituting  for  Y  =  ^32^21  w^ere 

k 


•32 


32 


A31  +  k31  +  k32 


1  - 


k  k 
*23  32 


(A31+k31+k32) {A21+k21+k23) 


(13) 


and 


21 


21 


A21  +  k21  +k23 


(14) 


gives 


20 

1-0 


4tt 


21  hV12  lV 


1  + 


A21+k23+k2.l 


"32 


8f  v 

1+  — [1+ 

83  \ 

J  13 


(13) 


X^2  Is  the  crossing  efficiency  for  populating  level  2  from  level  3  and  y is 
the  radiative  efficiency  for  20  luminescence. 


Case  IB(l-O).  With  the  same  assumptions  as  previously  given  (k^^ »  ^3* 
A32^^  *  an<^  assuming  radiative  excitation  of  level  2  with  non-radiative 
crossover  to  level  3  and  then  radiative  deexcitation,  the  ratios  of  ^/n^ 
and  are  &iven  by 


°1 


12 


E  +  E* 
v  V12  V12 


(16) 
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n3  c  K2 

ni  ”  8, 


1  t  V12  V12 


E.  +  E*  A„  +  k„,  +  k 


‘31  31  32 


where  all  terms  are  defined  in  the  manner  previously  given  for  similar  or 

* 

identical  terras.  The  term  E  is  defined  as 

V12 

*  -  CA21  _  82CA21  ^  ( 

V12  Y21E21  81Y21B12 

and  defined  as 


Y21  = 


A21  21  k23  ^A21+k21+k23) (A31+k31+k32) 


k23k32 


Case  I  Ba(l-»2,  2+1) .  Solving  for  in  terms  of  n^,  utilizing  equations 
(16)  and  (17)  and  substituting  for  n^,  gives 


,  (V 


A21  hv12  "T 


A31+k31+k32  g2 


+  —  1  + 


Case  I  Bb(l+2,  2+3,  3+1) .  Solving  for  n^  in  terms  of  rtj,  utilizing  equations 
(16)  and  (17)  and  substituting  for  =  X23?3i  where 


A21+k21+k23  1  ' 


_ k23k32 _ ‘ 

(A21+k2i+k23) ^A31+k3i+k32) 


3!  A31  +  k31  +  k32 


gives 


3F_  .  =  A31  hv13  nT 

3+1 


A31+k31+k32 


81  vi  2 

l  +-i  (l  + 

R2  Ev 

i  v12 


X23  and  y33  are  exactly  analagous  to  x32  anc^  ^21  resPect^v‘i^y ' 
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Case  IIaU--^).  In  this  case  non-radiative  rate  constants  and 
can  be  assumed  negligible  and  radiative  rate  constant  A^^O.  Assuming 
radiative  excitation  of  level  3  via  level  1,  the  ratios  n3/n^  and 
are  given  by 


_^3  =  £3 
nl  61 


'  Ev 


13 


E  +  E* 
V13  V13 


21 


13 


n„ 

- E* 

2 

+  —  k 

1+  v“ 

32  32 

83  12 

E 

V13 

is  defined  by 

r*  **31 

C83A31 

V13  Y31B31 

S1Y31B13 

defined  by 

A 

V 

31 

31  A3l  +  A32  +  k31  +  k32 


(24) 


(25) 


(26) 


(27) 


and  Y32  by 


'32 


32  A31  +  A32  +  k31  +  k32 


(28) 


11 


where  the  cyclic  crossover  terra  for  such  a  case  is  omitted,  i.e.,  the  upper 
and  middle  levels  are  not  both  radiatively  connected  to  the  lower  level 
(also,  k^3  and  k23  are  negligibly  small).  The  terms  k23  and  k^2  are  the  non- 
radiational  rate  constants  for  collisional  de-excitation  and  excitation  and 
are  related  by 

Si  -h* (VtT)  (29> 


where  E^2  in  the  energy  separation  between  levels  1  and  2,  T  is  the  flame 
temperature,  and  k  is  the  DolLr.mann  constant. 
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Imm-V «.  ... 


Case  II  Aa(l-»3,  3-*l) .  For  this  case,  solving  for  n_^  in  terms  of  n^., 
utilizing  equations  24  and  25  and  substituting  for  Y  from  equation  26 


gives 


47  A31  hv13  "t 


1  +  —  1  +  — 

8-3  E, 


i1  +  T±  e*p(-Ei2/kT>j  +  [ 


A32+k32 


Case  II  Ab(l-»3,  3  >-2) .  For  this  case,  solving  for  n^  in  terms  of 
utilizing  equations  24  and  25  and  substituting  for  Y^  gives 


BF„  „  =  4rr  A32  hV23  "t 


1  +  —  1+7 


—  1+  -=•  exp(-E12/kT)  + 

v13'  1  J 


A32+k32 


Case  II  B  (2-»3).  For  this  case,  non-radiative  rate  constants,  k^3  and  k23 
are  negligible  and  the  radiative  constant  A^^O.  Assuming  radiative  excitation 
of  level  3  via  level  2,  the  ratios  n^/n^  and  n^/n^  are 


n3  g2 


k21  exp(~E12/l:T) 


5,  E* 

3.  8  o  ^  o  o 

i"k21  f1+—l  +A31  +  k31 
3  V23 


^3  =  ^3  V23  (3- 

n2  "  g?  Ev  +  EC 

1  (.  v23  23  ' 

Case  II  Ba(2-*3,  3->2).  Solving  for  n^  in  terras  of  utilizing  equations 
32  and  33  and  substituting  for  Y.^  utilizing  equation  28  gives 


4 it  A32  hV23  "t 


8-,  v 

1  +  —  1  +  ~ 
8  a  E 


1  +  —  exp  (E  /kT) 
8?  12 


,  A31+k3ll 

k12  J 


Case  II  Hh  3-»] ) .  Solving  for  n^  in  terras  of  n^,  utilizing  equations  32 

and  33  and  substituting  Cor  Y  ^  utilizing  equation  27  gives 


5  5 


3+1 

2-*-3 


[/iir]  A31hV13nT 


1  +  — 

83 


E* 

v 


1  + 


23 


'23 


1  +  exp(E12/kT) 


Sl+k3 


I] 
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Three  Level  System — Limiting  Expressions 

In  this  section  of  the  manuscript,  the  limiting  expressions  for  the  8 
cases  above  will  be  given.  The  limiting  expressions  correspond  to  the  low 

: k 

intensity  case  (conventional  sources,  i.e.,  E^  >>  E^)  and  the  high  intensity 
case  (laser  sources,  i.e.,  E^  >>  E^  ). 


Case  I  Aa(l+3,  3-*-!) .  Low  Intensity. 


B„  (I.o ) 

F3+l 
1+3 


=  [t]  A31hV13nT  (^) 


\h) 

V13) 

UJ 

IE* 

L  v,-' 

13 


(12a) 


Case  I  Aa(l+3,  3+1) .  High  Intensity. 


B  (Hi) 
3+1 
1+3 


A31hV13nT 


1  +  —  + 


32 


S3  A21  +  k23  +  k21 


(12b) 


Case  I  Ab(l+3,  3— »2,  2+1) .  Low  Intensity 


B_  (Lo) 
F2+l 
l->3 


S3k32Ev..  I 

z 

A21hV12nT 

13 

(attJ 

g  (A  -fk  +k  )E* 

1  21  23  21  v23„ 

(15a) 


Case  I  Ab(l+3,  3 — *2,  2+1).  High  Intensity. 


V,0"  .  te 

l->-3 


A21hV12,1T 


1  + 


A21+k23+k2l) 


l  ^32  i 


1  +  — 

Si. 


(15b) 


C a se  I  Bn (1+2,  2+1).  Low  Intensity 


Bf2>1(Lo)  =  [/.it]  A2IhvI2nT 
1+2 


r  r’2Ev 


12 


'*1  V, 


(20a) 


(35) 
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3.  Tn  saturation  fluorescence  measurements  of  n^,  Bp  under  steady 
state  conditions  must  be  measured,  i.e.,  the  radiative  lifetime  of  the 
excited  state  must  be  much  less  than  the  temporal  pulse  width  of  the 
exciting  source.  Also,  the  source  intensity  must  be  such  that  the  fluores¬ 
cence  radiance  will  increase  by  less  than  a  few  percent  if  the  source 
intensity  increases  by  orders  of  magnitude.  Pulsed  resonance  spectroscopy 
where  saturation  is  achieved  and  the  transient  decay  is  observed  to  determine 
n^  is  difficult  to  perform  experimentally. 

4.  For  3-level  Case  I  atoms  and  for  1-3  excitation,  Y  contains  a 

r 

cyclic  term  involving  levels  2  and  3,  i.e. ,  crossover  and  crossback. 

This  term  is  identical  to  the  intersystem  crossing-recrossing  tern  in 
phosphorescence  of  organic  molecules. 

5.  Spatial  atom  profiles  (determination  of  absolute  n^  values) 

via  saturation  fluorescence  spectrometry  (laser  excitation  with  observation 

of  fluorescence  from  a  small  volume  excited  by  the  laser)  can  be  obtained 

if  the  atom  (or  molecule)  is  excited  and  de-excited  as  a  2-level  system, 

or  if  the  atom  (or  molecule)  is  a  3-level  Case  I  (levels  2  and  3  close 

together).  In  the  former  case,  one  only  has  to  measure  Bp  and  if 

V12’  ®1  an<^  3re  ^ien  can  t>e  calculated.  In  the  latter  case, 

one  has  to  measure  (absolute  radiance)  3-1  fluorescence  3—1  and  the 

Bp  Bp  1-3 

ratio  n./n_  by  measuring  the  ratio  2-1/  3-1  [=  n„/n.  =  ^-/(A-.+k-.+k,,..)  ] . 

j  l  1_^3  X-v3 

By  use  of  equation  12b,  assuming  saturation  (high  intensity  source)  has 
been  achieved,  n^  can  be  estimated  by  substituting  of  values  for  A^, 
hVf^,  gp,  and  the  ratio  n2/n3>  For  atoms  of  type  Case  II,  it  is  not 
possible  to  determine  n_  unless 
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14 


•i  -y-  ,/iTj 


resulting  in  the  simplification  of  equation  30b  to 


ii  (Hi)  ='= 
3->l 
1-+3 


A.,  ,  bv 


f. 


31  13  T  [nl  +  fij 


or  simply  the  2  level  approximation  for  a  3-level  case. 

6.  The  definition  of  n  saturation  spectral  i r rad la nee  is  possible, 
i.e.,  the  source  spectral  irradiance  chore  the  fluorescence  radiance  becor.cs 
503  of  the  maximum  possible  value  but  no  meaningful  information  results  for 
the  3-level  system  since  E  depends  upon  not  only  measureahle  or  known 
radiative  transition  probabilities ,  A' s , hut  also  upon  difficult  to  measure 
and/or  unknown  radiationless  rate  constants,  k's.  Therefore,  in  this  note, 

no  attempt  was  made  to  give  such  values  despite  the  simplicity  of  the  resulting 

B  expressions.  (Refer  to  Appendix  B  for  a  general  approach  to  evaluation  of  Es) 
r  V 

7.  The  relative  error  resulting  in  the  evaluation  of  if  the  atom 

is  assumed  to  be  a  2-level  case'*'  and  it  is  in  fact,  a  3-level  Case  I  (Na)  is 

-45%  (excitation  of  level  2)  or  -22%  (excitation  of  level  3)  and  if  the  atom 

is  assumed  to  be  a  3  level  Case  1  hut  the  crossover  terms  involving  3-<-+-2 

are  neglected^  is  :  +44%  (excitation  of  level  2)  and  +11%  (excitation  of  level 

2 

3).  These  values  are  calculated  for  an  air-acety lore  flame. 

8.  Apparently  anomolous  saturation  spectral  irradinneo  values  will  he 
obtained  if  the  source  irradinneo  is  not  constant  across  the  focal  volume. 

Daily*'’  has  discussed  the  effect  of  a  Gaussian  laser  beam  ivradiance  dis¬ 
tribution  on  atonic  fluorescence  signals  for  several  excitation  geometries. 
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APPENDIX  A1 


STEADY  STATE  APPROACH  TO  CONCENTRATION  EXPRESSIONS 

Assuming  that  the  rate  equations  approach  can  be  safely  applied  and 
that  the  3  levels  can  all  be  coupled  radiativeiy,  the  general  expressions 
for  the  rate  of  change  of  concentrations  in  states  1,2,  and  3  are  given  by 


.  B_„E 

dt  c  c 


+  A31+A32+k31+k32 
Eqn.  (Al) 


n3  + 


B_  ~E 
13 

- —  +  k., 

c  13 


nl  + 


B  _E 
23  V23 

- —  +  k,, 

c  23 


^2 

dt 


B23Ev„,  B21Ev 


23 


c  12  +  A21+k21+k23 


n2  + 


B32Ev 

“^-32-32 


n3  +  l 


Eqn.  (A2) 

"T  =  nl  +  n2  +  n3 


fB12E,i2 

— - —  +  k._ 
c  12 


(A3) 


By  setting  the  above  equations  equal  to  zero  for  steady  state  conditions,  we 
can  solve  for  the  ratios  n2/n3,  and  n^/n^  obtaining  the  following 

expressions 


( B,  E 
12  v 


12 


+  k 


12 
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B,  E  B,„E 
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Obviously  (see  specific  cases  in  the  manuscript)  these  expressions  simplify 
when  one  particular  3-level  case  is  considered  and/or  when  excitation  is 
performed  at  a  selected  transition. 

For  an  atom  with  a  complex  set  of  energy  levels  and  potential  transitions, 
e.g.,  the  transition  elements,  the  actinide  elements,  and  the  rare  earths, 
even  more  complex  expressions  result  because  more  than  3  levels  will  be 
involved.  In  such  a  case,  summations  over  activation  and  deactivation  routes 
from  or  to  a  specific  level  must  be  made.  Because  of  the  limited  use  of  such 
structurally-spectrally  complex  atoms  in  flame  diagnostics  via  atomic  fluores¬ 
cence,  no  specialised  expressions  will  be  given  here. 
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APPENDIX  A2 


SATURATION  SPECTRAL  IRRAD LANCE,  E® 

*  v 

2—5  8  9 

Omenetto  et  al.,  D'Olivares,  and  Daily  have  utilized  a  parameter 
called  the  saturation  spectral  irradiance,  E^,  which  is  useful  experimen¬ 
tally  because  it  is  the  value  of  source  spectral  irradiance  when  the 
fluorescence  radiance  is  50%  of  the  maximum  possible  value  (i.e.,  B_  is 

w 

plateau  value  when  E  — *») . 

v 

g 

For  a  2-level  atom,  E  is  given  by 

V12 
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where  the  terms  are  as  defined. previously  (i.  is  the  ground  state,  u  is  the 

radiationally  excited  upper  state,  and  u*  is  the  other  excited  upper  state). 

Experimentally  one  can  plot  Bp  vs  Ey  and  determine  E^  by  simply  estimating 

E®  at  By  12.  Theoretically,  it  is  difficult  (in  many  cases,  impossible) 
max 

to  estimate  Es  because  of  the  lack  of  information  on  the  radiationless  rate 

constants,  ks.  The  expression  for  Es  converts  to  the  E®  expression  for 

V£u  12 

the  2  level  atom  if  one  assumes  all  k  and  A  terms  with  u'  are  zero. 

For  a  three-level  atom  of  the  t£  type,  a  different  expression  results 


where  u  is  the  upper  level,  l  is  the  lower  level  from  which  radiational 
excitation  occurs,  and  £'  is  the  other  low  level  (for  TZ,  u  =  3,  l  *  1  or  2, 
and  l'  =  2  or  1,  respectively).  The  above  expression  (B3)  converts  to 
the  two-level  case  with  appropriate  substitutions. 
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ABSTRACT 


Five  methods  for  measurement  of  flame  temperatures  based  upon 
laser  excited  fluorescence  are  discussed  with  suspect  to  the  derived  ex¬ 
pressions,  the  assumptions  necessary,  the  advantages,  and  the  disadvant¬ 
ages.  The  use  of  each  of  the  five  methods  for  measuring  spatial  (flame 
volumes  of  1  mm3)  and  temporal  (flame  temperatures  corresponding  to 
1  ys  time  period)  is  given.  The  five  methods  consist  of  one  based  upon 
linearity  between  the  fluorescence  signal  and  the  laser  spectral  irrad- 
iance ,  of  three  based  upon  saturation,  and  of  a  fifth  which  is  not 
critically  dependent  upon  the  laser  spectral  irradiance. 
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I  nt  roduc  t ion . 


Atomic  fluorescence  excited  with  a  spectrally  quas i -cont inuum  source 
has  been  shown  to  be  a  reliable  method  of  measuring  flame  temperatures.*  2 
In  fact,  it  combines  the  advantages  of  optical  methods,  namely  not  disturb¬ 
ing  the  combustion  process,  with  the  capability  of  providing  spatially 
resolved  measurements:  the  latter  feature  is  common  to  all  "scattering" 

O 

methods,  including  those  based  on  Raman  methods. 

The  use  of  a  laser  as  an  excitation  source  for  atomic  fluorescence 
has  the  obvious  advantage  that  the  spatial  resolution  can  be  improved  be¬ 
cause  of  the  low  divergence  of  the  collimated  beam  coupled  with  its 

high  spectral  radiance.  The  use  of  a  pulsed  laser  has  the  additional 
characteristic  when  used  with  a  gated  detection  system  of  measuring  tempor¬ 
al  (short  time  period)  temperatures. 

9  10  11 

Both  Measures  and  Omenetto  and  Winefordner  ’  have  discussed  the 

use  of  laser  saturation  of  probe  transitions  and  the  measurement  of  appro¬ 
priate  fluorescence  transitions  for  estimation  of  flame  temperatures. 

These  reports  indicated  that  neither  calibration  of  the  electro-optical 
detection  system  was  needed  nor  was  measurement  of  the  laser  power. 

In  this  paper,  we  give  the  theoretical  expressions  of  5  fluorescence- 
based  methods  for  flame  temperature  measurement  including  predicted  ad¬ 
vantages  and  disadvantages  of  each  method.  Two  of  the  methods  have  been 
previously  described.*  **  All  five  methods  are  described  with  consisten 
terminology  using  the  well-known  steady  state  approach;  the  general 
expressions  and  symbolism  are  according  to  the  manuscript  of  Boutilier, 
et  al.*2 
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Theoretical  Basis 


Boundary  Conditions  and  General  Hxprossions.  In  all  five  methods, 
a  3-level  probe,  such  as  T1 ,  In,  Ga ,  Pb ,  etc.,  is  assumed  to  be  in  the 
flame.  This  can  be  experimentally  accomplished  by  nebulization  of  a 
salt  solution  containing  the  inorganic  probe  to  produce  an  aerosol  which 
is  partially  and  reproducibly  transferred  to  the  flame  system  under  study. 
Such  methods  are  described  by  Herrmann  and  Alkemade.^  It  can  also  be 
accomplished  by  metal  vapor  generation  or  vaporization  of  volatile 
metal  complexes. 

An  energy  level  diagram  with  appropriate  symbols  for  activation  and 
deactivation  rates  is  given  in  Figure  1.  For  an  atomic  probe,  such  as 
Tl,  excited  at  either  l->-3  or  at  2> 3  (see  Figure  1)  the  following  assump¬ 
tions  can  be  made:  %  0 ’  k23  %  0;  A21  =  0.  Solving  for  the  ratio  of 

n,/nj  and  f°r  steady  state  conditions  “  (dn^/dt  =  0  and  dn2/dt  =  0) 

under  continuum  or  pseudo  continuum  excitation  gives 


The  fluorescence 


Using  equations  (1)  and  (2)  with  the  necessary  condition 


nT  =  nl  +  n2  +  n3  W 


and  with  either  Equation  (4)  or  (5)  allows  evaluation  of  BP  and  Bp 

r3-l  r3-2 


in  terms  of  n.j.  and  appropriate  excitation  conditions  i.e.,  Ey. 


1  -  7 

Method  1.  Two  Line  Linear  Method?  In  this  method,  the  fluorescence 

at  3-1  excited  via  2  11 3 ,  B,.  ,  and  the  fluorescence  at  3-*- 2  excited  via 

h3-l 

2-3 

1-3,  Bp  are  measured  under  1  inear  Bp  vs  conditions.  Evaluating 
o-l 
1-3 
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liquations  4  and  5  for  these  conditions  gives 


*-7tT-*  A31hv31nT 
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>,E*  1  +  -1 

“  v23  82 


exp  C E  j 2/kT) 


^ 4tt ^  A52hv32nT 


'  C1  *  - 

'i-L  e  l 


exp (E , 9/kT) 


where  all  terms  are  defined  in  the  caption  of  Figure  1  except 


E  =  modified  saturation  spectral  irradiance  with  excitation  at  v 


c  A32  87rhv23 

H,,  =  „ — v —  =  — t — • —  and  Y_ 


v 2 5  B32A  32  c2Y.7  *2  A31  +  A32  *  k32  +  k32 


=  modified  saturation  spectral  irradiance  with  excitation  at  v 

v  1  3 


*  c  A~^  Snhvt^  A-^ 

1 3  B 3 1 Y 3 1  c2Y5]  3nd  Y'^  A3 1  +  A32  +  k31  +  k32 


where  v  is  the  frequency,  in  lie,  of  the  appropriate  excitation  transition 
and  Y  is  the  quantum  efficiency  of  the  designated  fluorescence  process  of 

the  probe  in  the  specific  environment  of  concern.  Dividing  Equation 

- 1  -  1  -  l 

(8)  by  (7)  and  using  t he  more  common  units  J  s  in  “nm  tor  the  source 
spectral  irradiances,  we  obtain  the  following  expression  for  the  temperature, 
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l:-r,A 

= 

i;x,-  X,, 

?)U|r-=^  )  +  6  ?.n(^-)  +  g-n  -■— ) 

'A  ,  -  13  F,  , 


Method  II.  Saturation  Two  Level  Peak  Method  (Pseudo-Two  Level  Steady  State) 

In  this  method,  the  ratio  of  the  fluorescence  radiance  resulting 

from  pumping  both  1 '5  and  2 -*-3  sequentially  is  measured.  However,  in 

t  it  i s  method,  it  is  assumed  that  at  a  sufficiently  high  source  spectral 

irradiance  (and  fast  rise  time),  the  pumping  time,  t^ ,  is  sufficiently 

1 

short  to  assure  that  t  >>  k~,  +  k,7  +  A-,  +  A-,.  In  this  case,  the 

p  ol  32  aloe 

system  reduces  to  a  pseudo-two  level  system  consisting  of  only  the  states 

B  i  a  F  ]  7  ^21  ^  1 

being  pumped.  In  other  words,  — 11  and  — dominate  over  all  other 

rates  in  a  period  of  time  which  is  much  shorter  than  the  upper  level's 
lifetime,  t.  In  this  case  not  only  must  t^  be  sufficiently  small,  but  the 
measurement  system  for  the  resultant  fluorescence  radiance  must  also  be 
sufficiently  fast  to  resolve  temporally  the  fluorescence  temporal  waveform 
in  order  to  obtain  accurate  measure  of  the  steady  state  levels  for  the 
fluorescence  radiance,  which  first  rises  to  a  pseudo-two  level  steady-state 
value  (peak)  with  a  response  time  controlled  bv  t  and  then  relaxes  to  its 

p 

three  level  steady  state  value  controlled  by  k-j,  k~7,  k]7,  k-,^,  A-^,  A_, 
(see  Figure  2).  This  process  has  been  di scribed  by  Measures",  and  by 
Omenetto  and  U’inefordner .  '  ^  ' 


f  The  pumping  time,  t^,  is  given  by1  ;  t  =  t 


+  1 ; " 

v0  VC 


—I  ,  where  t  is  the 


effective  lifetime  (t*:  1")  and  F.?  :%  For  C  >>  HS  ,  then  t  f  f,',(Y). 

sp  V()  \  v0  v0*  P  nV 


*?  O 


The  3*1  fluorescence  radiance  for  1*3  excitation  and  for  2*3 
excitation  under  saturation  2-level  conditions  are  given^’^  by 


BF3^  A31hv31n3  A31hv31nTfg1  "+  g'”^  (12) 

1*5 


3*1 

2*3 


A 


hv-T  ru 
31  ol  :> 


^ 4 7r ^  A31hv31nT^g1')  exP(-Ei2/kT)  ( 


§3 


g2  +  §3 


where  n~  is  essentially  the  population  of  the  ground  level  prior  to 
laser  excitation 

Dividing  Equation  (12)  by  (13)  gives 


(13) 


3*1 

1*3 


=  ( 


e2 


+  g  ~ 


3*1 

2*3 


— )  ( — )  exp (E, , 
83  ^2  12 


/kT) 


(14) 


Solving  for  T  and  introducing  the  g  values  for  T1  gives 


E^/k 


‘II  ■ 

nwf 

£n 

4 

1*3 1 

3 

bf 

3*i  / 

mi 

i  2*3/ j 

(15) 


Method  III .  Saturation  Two  Line  Method  With  Sequential  Pumping . 

Tn  this  method,  the  ratio  of  fluorescence  signals  are  measured  after 
sequential  excitation  of  fluorescence  under  saturation  conditions  assuming 
the  system  has  relaxed  to  its  three  level  state  (see  Method  II).  Because 
several  ratios  can  be  used,  four  fluorescence  radiance  expressions  will  be 
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given  by  evaluating  Equations  4  and  S  assuming  saturation  conditions 
apply,  i.e., 


S-i ' (i7)  A^h”3inT  ir*-*r~7r>  ,  „  ~  .  *3y*v2 

1-3  L  ~ gj  *  (17>  exp(-Ei2/kT)  *  k21-  ■ 


32  32 


v31A31  J  F3-l 


g  =  (r~)  A,1h’j71nT  - - - - - - - - r —  (IS) 

F  W  31  31  T  g2  -  g3  ,t.  ~  _  ,  a31  +  k31 

2-3  ~~g~z  +  (g^}  eXp(t!2/kr)  +  ~ T ~J 


Bp  =  hrV1  BF 
b3-2  ,  V31A3  1  /  h3-l 


In  general,  the  ratio  of  Bc  /Bc  ,  where  i  and  j  correspond  to  levels 

h3-l  1  3- j 
1-3  2-3 

1  or  2 ,  is  given  by 


\  g2  +  S3  Si  8l,A31k31 J 

iSd  -iT"  ‘  h  '  exp(^ 


v, .A,  . 
■5J  3j 


g3  A32  +  k32 


as  long  as  a  probe  is  chosen  where  cxp( -E j -,/kT)  <<  1  (or  E^,/ T  <<  10  cV  K  ) 
this  corresponds  T1  or  Pb  for  flame  temperatures  below  'vlOOO  K  but  does 
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no t  apply  to  (la  or  In  for  llano  temperatures  between  2000  and  4000  K. 
Solving  for  T  (and  introducing  g ^  =  g-  =  1  ;  g2  =  2)  gives 


El2/k 


III 


£nl 


?  + 


j2  j2 


r 


T- 


21 


A, ,  +  k- 
al 


r 


aj  ad 


a+ii 

l-3( 


b,w 

21  ^  71  2-3 


^31  +  ^31 

2  + 

k21  -J 

(21) 


If  in  addition,  the  following  assumption  is  valid,  k21  >>  (A.^  +  k31)  ^ 
(A--,  +  ^3?)  >  then  Equation  (21)  reduces  to 


El2/k 


III 


in 


-  t. 

2(^3J^)('-Tri^) 


1  v  ,  .  A  .  . 

3  i  a  l 


-  3 


3-h 

2-3 


(21  ) 


Finally  if  n^,^  30  n-j.  and  since  v3^A3j^v3i^31  %  *  >  then  Equation  (18) 
reduces  further  to 

(21"  j 


Method  IV.  Saturation  Two  Line  Method  With  Simultaneous  Pumping. 

In  this  method,  which  is  similar  to  Method  III,  the  ratio  of  the 
fluorescence  radiance  resulting  from  pumping  both  1-3  and  2-3  simultaneously 
to  the  fluorescence  resulting  from  pumping  one  of  the  fluorescence  transitions 


TIII  ”  % 


E12/k 


£n 


/B 


F3-i^ 

1-3 


3-1 
2-3 
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given  in  Method  II  (under  saturation  conditions)  is  measured.  From 
liquations  (1),  (2),  (4),  (5),  and  (6),  the  fluorescence  radiance  for 
simultaneous  pumping  is  given  by 


BF-  =  hv 

O-l 

<^3> 

2-3 


83 


31  31  T|  g,  ♦  g2  ♦  g3 


(22) 


if  the  fluorescence  at  3-1  is  measured  and  by 


1  3-2 
2-3 


v  ,  A  ~  n 

bf 

v31A31 


3-1 


(23) 


<\:i> 


if  the  fluorescence  at  3-2  is  measured.  For  the  case,  where  the  ratio  of 

Bp  to  Bp  is  measured 
*'5-1  j-1 

i  .  r  nr 


3-1 

<\+\> 

2-3 


g2 


+  «3 


«1 

83 


.  gl  /31  +  k3 1 

h  (~  "li 


exp(E12/kT) 


BF 


3-1 

2-3 


8l  + 


82  +  83 
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(24) 


Solving  for  T  from  Equation  (24)  and  evaluating  the  g's  gives 
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and,  as  previously  in  Method  HI,  if  k 


and  finally  as  in  Method  III,  if  n2''<  30 


then 


(25 


The  case  where  the  ratio  of  B 


lp  to  Bp  is  measured  is  not  considered 
^  3->l  ^3-1 

<&> 


here  because  this  case  is  limited  by  the  precision  in  the  measurement 

c.  / .  c 

fi.e.,  the  ratio  varies  from  Wd)  at  T-»-0  to  a  maximum  of  S^T)  at  T  -»•  «)  . 


In  the  above  cases  in  this  method,  it  is  assumed  that  measurement  of 
the  fluorescence  radiance  occurs  in  the  three  level  temporal  regime  (see 


Method  II).  If  instead  we  take  the  case  where  the  2  3  excitation  occurs 


just  prior  to  the  simultaneous  pumping  of  1  3  and  2  3  (see  Figure  3)  and 


if  t  is  fast  enough  to  achieve  a  pseudo-two  level  system,  then  the 

P  >  s 


observed  fluorescence  radiance  ration  B~eak/B 

3-*-l 
2-*3 


r  is  given  by, 

3-*l 

l-*3 

_> 


3-*-l 

<;■*?> 


Si(g2  +  g3) 


— -  - ±— - r  exp(-E. ,/kT) 

peak  g2(gi  +  g2  +  g3 ))  12 


(26) 


F  -  , 
j->-1 

2->-5 


Substituting  for  the  g's  yields 


E,  ,/k 


IV 


peak 


5,n 


r 

/bf  \1 

/  3->l  \ 

<1-3 

<!> 

<2-5> 

Vl  / 

. 

Hi 

(27) 


It  should  be  noted  that  Tjy  is  completely  indipendent  of  k^2,  ^21’  k5l » 

peak 

k32 ’  A 5 1 ’  A32- 

Method  V.  Fluorescence- thermally  Assisted  Fluorescence  Method. 

In  this  method,  only  one  laser  excitation  wavelength  is  needed  and- 
achievement  of  saturation  i_s  not  necessary  (see  Figure  4).  Here,  the 
upper  level  of  the  probe,  3  in  our  case,  is  excited  via  a  laser  and  then 
either  the  ratio  of  the  fluorescence  3->-l  to  the  thermally  assisted  fluores¬ 
cence  from  a  close  lying  upper  level,  say  4  +  1,  is  measured  or  the  ratio  of 
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the  thermally  assisted  fluorescence  from  two  close  lying  upper  levels, 
say  S-*l  to  4-*-l  is  measured. 

Solving  the  steady  state  equations1^1  (in  a  system  similar  to  that 
shown  in  Figure  4)  for,  n^/n^,  gives 


^3  =  k43  +  k41  /  A4 1 
n4  ^34  k34 


(28) 


where  it  is  assumed  that  A^~  =  0  and  that  the  color  temperature  of  the 

sources  of  excitation  are  large  compared  to  the  thermal  temperature  under 

consideration;  this  assumption  is  valid  under  most  cases  employing  laser 

B 1 3  E  i  5 

excitation  (i.e.,  - — -  >>  k ,  ~ ,  k14,  or  k34).  The  ratio  of  fluorescence 

radiances  Bp  /BTp  is  proportional  to  n-/n^.  In  addition  if  we 
assume  that  the  energy  difference  between  levels  4  and  3,  E-4 ,  compared 
to  Eu  ,  is  significantly  small  (i.e.,  exp(-F-34/kT)  >>  exp(-E15/kT) 

to  allow  us  to  assume  that  k^4  >>  k41  +  A41,  then  equation  (26)  reduces 
to, 


k43 
"1  *34 


(29) 


which  is  simply  the  Boltzmann  distribution  where. 


^4  ®4  "^34^kT 

n3  =  «3  C 


(30) 


Including  other  levels  such  as  in  the  Tl-type  atomic  system  will 
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further  complicate  analysis  of  n^/n^  f°r  various  levels,  but  under  the 
same  set  of  limiting  conditins,  where  it  is  assumed  that  the  Bolztmann 
distribution  applies  (i.c.,  Ei  not  t00  large  for  the  given  T)  then  wg 

j 

will  employ 


\  .  «i  e‘Eji/kT 

11  i 


(31) 


for  all  levels  with  i  and  j a  the  upper  level  being  pumped. 

In  the  case  of  Tl-like  atomic  probes  we  can  obtain  the  thermally 
assisted  fluorescence  radiance  expression  given  by: 


1-F^i  =  ^  Ailhvil(^  exp(-Ei3/kT)nJ 


i-l 

1-3 


gl  +  g3  81 
___  +  _  exp(-E12/kT)  * 


A32  +  k32 


(32 


21 


g; 


RrC.  .  =  (£-)  A.-h  )  exp(  E.,/kT)ru.  - - - - - 

TFi-1  W  ll  ll'-g,'  i3  ng.tg,  g. 

2-3  -4 - -  +  ^  exp(-E12/kT)  + 


g3 


g3 


12 


(35 


Taking  the  ratio  of  Bn  to  BTC  ,  where  h  is  the  lower  level  of  the 

^3-h  lt-i-l 

g-3  g-3 

fluorescence  transition,  g  is  the  lower  level  (1  or  2)  of  the  excitation 
transition,  and  i  is  the  upper  level  produced  by  thermal  excitation  of  the 
laser  excited  level,  3,  gives 
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'hi'oret  ion  1  Comparison  of  Methods  . 

let  hod  I .  The  advantages  of  this  method  are: 

(i)  saturation  is  not  necessary  and  in  fact  must  not  be  approached; 
thus  relatively  low  laser  irradiances  can  and  must  be  used; 

(ii)  it  is  independent  of  non- radiative  activation  and  deactivation 
rate  constants  and  radiative  deactivation  rate  constants,  i.e., 
no  assumptions  are  needed  to  obtain  the  expression  for  Tj  ;  and 

(iii)  temporal  flame  temperatures  can  be  obtained  by  passing  both 
beams  through  the  flame  of  interest  but  delaying  one  beam 
(pulse)  and  measuring  the  ratio  of  Bp  values  for  the  individual 
pulses.  The  time  integrated  as  well  as  the  time  dependent 
fluorescence  radiance  may  be  measured  in  this  method,  since  it 

c 

is  assumed  that  B.  E.  >>  k.  and  E  <<  E  (i.e., 
lu  lu  lu  Vq  v0 

c 

kul's  +  Aul ' s  >>  B1uE1u  ).  Therefore,  the  fluorescence 

- 

temporal  waveform  will  follow  the  source  temporal  waveform 
linearly  at  all  points  in  time  where  these  conditions  hold. 

Both  the  fluorescence  radiance  and  the  source  spectral  irradiance 
require  the  same  integration  time,  t  .  and  the  same 
time  delay,  xdelay 

The  disadvantages  of  this  method  are: 

(i)  because  the  fluorescence  radiance  values,  Bp,  depend  linearly 
upon  the  source  spectral  irradiance,  E^ ,  and  the  flame 
quantum  efficiency,  Y,  the  fluorescence  signal  strengths  are 
lower  than  in  saturation  based  methods  for  equal  n^'s; 

(ii)  for  the  3-»-l  fluorescence  process,  the  post  filter  effect 
causes  a  reduction  in  the  fluorescence  signal; 
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(iii)  Bp  vs  r-:  must  he  linear,  i.e.,  measurements  must  be  taken  on  the 

linear  portion  of  the  Bp  vs  plot; 

(iv)  the  source  spectral  irradiances  at  1+3  and  at  2-*- 3  must  be  cal¬ 
ibrated  in  relative  irradiance  units  and  the  fluorescence 
radiances  at  3+2  and  at  3+1  must  be  calibrated  in  relative 

radiance  units  assuming  equal  r  and  r,  ,  : 

M  meas  delay 

(v)  the  two  laser  beams  resulting  in  excitation  at  1+3  and  at  2+3 
must  be  matched  in  terms  of  optical  alignment  and  must  not 
cause  saturation  even  in  the  central  portion  of  the  beams; 

( v i )  single  pulse  temperature  measurements  require  an  accurate 
repeatable  measure  of  A^aser  from  pulse  to  pulse  since 

the  source  spectral  irradiance  must  be  accurately  known;  and 

(vii)  the  paired  measurement  of  1+3  and  2+3  excitation  must  occur 
on  a  time  scale  such  that  AY/ At  =  0. 

Methods  II  and  IV  eaj.-  The  advantages  of  these  methods  are: 

(i)  the  non- radi at ional  and  radiational  rate  constants  do  not 
enter  the  Tjj  relationship,  i.e.,  no  approximations  are 
requ i red ; 

(ii)  temporal  temperatures  are  more  easily  obtained  since  only  the 
fluorescence  radiances  need  to  be  measured;  and 

(iii)  high  speed  "real-time"  temperature  computation  is  more  accessible 
due  to  the  simplistic  relationship  for  T^j,  and  Tjy  aj_. 

The  disadvantages  of  these  methods  are: 

(i)  the  laser  pulses  must  he  spatially  homogeneous  which  usually 
requires  spatial'  filtering  of  the  dye  lasers  output  limiting 
tiiis  technique  to  high  power  Nd-YAG  and  excimer  pumped  dye  laser 
systems  for  low  quantum  efficiency  combustors  (atmospheric 


hydrocarbon  flames); 


(ii)  the  laser  pulse  rise  time(s)  must  be  sufficiently  "fast" 
enough  to  allow  t^1  to  be  much  larger  than  all  the  other 
collisional  rates,  which  may  further  limit  the  technique  to 
cavity-dumped  systems  when  Y  is  very  small; 

(iii)  the  detection  system  must  be"fast"  enough  not  to  significantly 
distort  the  temporal  information  content,  i.e.,  the  pseudo- 
two  level  peak  region,  and  have  sufficient  sensitivity  to 

7 

assure  a  single  pulse  S/N  >  10  (i.e.,  greater  than  10  photons 
measured  per  t„„„) ; 

(iv)  scatter  errors  (Mie  scatter  primarily)  may  be  appreciable  for 

resonance  processes;  non- resonance  fluorescence  minimizes 

this  problem,  it  is  recommended  that  Bp  and  BR  be 

1  3+1  ^3+1 

1  +  3  2+3 

used  in  Method  II  (and  III)  because  the  weaker  signal,  BP 

h3+l 

2+3 


is  immune  to  scatter  and  B 


c  (as  well  as  BP  i 
o+l  *'3+1 

1+3  <l+3> 

2+3 


in  Method  IV) 


is  (are)  large  compared  to  the  scatter  signals; 

(v)  the  post  filter  factor^’^  will  cause  errors  for  any  fluores¬ 
cence  process  involving  the  ground  state,  e.g.,  BP  or  BP 

*"3+1  r  3+1 
1+3  2+3 


(or  Bp  );  however,  by  measuring  the  fluorescence  at  the  same 
o+l 
1+3 
2+3 

wavelength  for  both  fluorescence  processes  in  the  ratios 
(in  Methods  II,  III,  or  IV),  the  post  filter  effect  cancels; 


(vi)  the  laser  beams  must  not  only  be  homogeneous  as  in  item  (i) 
above  but  also  must  ha  ■  the  same  cross-sectional  areas  and 
must  be  optically  aligned  to  interact  with  the  same  region 
of  the  flame  gases,  and  in  fact,  the  same  species  when 
making  temporal  temperature  measurements;  and 

(vii)  full  saturation,  i.e.,  E  >  10ES  ,  must  be  achieved,  i.e., 

v0  v0 

it  is  not  desirable  to  plot  Bp  vs  and  to  extrapolate  to 
the  plateau. 

Methods  III  and  IV.  The  advantages  of  these  two  similar  methods  are 

(i)  since  "saturation"  must  occur,  only  fluorescence  radiances  need 
to  be  measured; 

(ii)  "real-time"  temperature  computations  are  easily  accessable 
assuming  the  validity  of  the  Tjjj"  and  Tjy'  equations; 

(iii)  single  pulse  temperature  measurements  are  easily  achieved  in  both 

methods;  however,  Method  IV  is  more  easily  adapted  since  the 

first  laser  can  be  set  at  E  and  Bn  measured  and  while 

v25  3-1 

2 -+3 

the  E  pulse  is  still  on,  the  second  laser  at  E  can  be 
v25  v13 

pulsed  yielding  Br  ,  i.e.,  there  is  no  need  for  delaying 

3-»  1 

<]^3> 

2-3 

one  pulse  relative  to  the  other  pulse  long  enough  to  be  certain 
that,  the  system  under  study  has  regained  thermodynamic  equili¬ 
brium,  that  the  fluorescence  temporal  waveforms  do  not  overlap 
and  that  all  electronic  echos  and/or  reflections  from  the  first 
pmse  have  been  damped  out  as  is  necessary  for  Methods  I,  II, 
and  !  t  I  . 
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The  disadvantages  of  these  methods  are: 

(i)  the  same  ones  as  (i),  (iv),  (v) ,  (vi)  of  Method  II; 

(ii)  the  laser  pulse  must  be  temporally  long  enough  to  assure  the 
pseudo-three  level  steady  state  condition; 

(iii)  while  extrapolation  to  the  maximum  fluorescence 

radiance  may  be  employed  in  Method  III,  it  can  not  be  used  in 
Method  IV;  in  addition,  it  can  only  be  employed  in  those  cases 
of  Method  III  where  steady  state  is  achieved  (i.e., 

^ laser  >:>  ZXY^ — +~~K — f  or  e*se  extrapolation  will  yield 

an  erroneous  result; 

(iv)  if  the  laser  pulse  and  measurement  times  become  too  long,  ther 
excited  state  reactions  must  also  be  taken  into  account,  and 
the  rate  equations  will  require  subsequent  modifications  (it 
should  be  noted  that  unless  photoionization  occurs,  Methods 

II  and  I V p c a p  do  not  have  this  problem); 

(v)  the  disadvantage  listed  in  item  (iv)  indicates  limited  use 

of  high  power  C.W.  or  long  pulse  flashlamp  pumped  dye  lasers 

while  the  disadvantages  in  items  (ii)  and  (iii)  indicate 
the 

limited  use  of  /^shorter  pulse  N7,  Nd-YAG,  and  excimer  pumped 
dye  lasers; 

(vi)  the  measurement  system  need  not  be  as  "fast"  as  is  required  i 
Method  II  but  it  still  must  be  able  to  measure  the  three  leve 
fluorescence  without  temporal  distortion  and; 

(vii)  unless  the  assumptions  made  to  obtain  Tjjj"  or  ^ni"  ant^  'l\ 
or  Tjy"  arc  valid  or  alternately  verified  then,  the  non- 
radintional  and  radiational  rate  constants  (i.e,  k's  and  A's) 
must  tie  known  to  obtain  an  accurate  temperature  (the  assumpt 
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made  to  obtain  T^,  and  T^j,  appear  to  be  reasonable  based 
upon  measured  values  of  k's  and  A's.^ 

Method  V.  The  advantages  of  this  method  are: 

(i)  saturation  is  not  necessary  but  can  occur; 

(ii)  there  is  only  one  laser  beam  (wavelength)  and  so  there  are  no 
matching  of  beam  problems  and  homogeneity  of  the  beam  is  not 
critical ; 

(iii)  scatter  and  post  filter  problems  are  minimal  because  the 

1+3 

value  is  large  and  immune  to  the  post  filter  effect  and  scatter 

errors,  whereas  BTC  is  immune  to  scatter  and  B-,,  is  also 
TFi+l  TFi+2 

fairly  insensitive  to  the  post  filter  effect; 

(iv)  temporal  temperatures  are  most  ideally  done  with  one  laser 
pulse  since  no  timing  problems  are  necessary,  but  rather 
simply  the  simultaneous  measurement  of  the  fluorescence  and 
emission  or  of  two  emission  signals; 

(v)  this  method  also  adapts  itself  readily  to  real  time  since  only 
the  fluorescence  radiances  need  to  be  measured  and  Ty  allows 
easy  temperature  computation  from  single  pulse  measurements;  and 

(vi)  as  in  Method  I  (for  Ey  J  ),  either  integrated  or  time 
dependent  measurements  are  applicable. 

The  disadvantages  of  this  method  are: 

(i)  the  collisional  rate  constants  k~.  and  k~,  must  be  sufficiently 

OX  O  J 

large  to  assure  that  levels  i  and  j  are  thermally  populated  on  a 
time  scale  of  that  of  the  laser  pulse;  temporal  analysis  of 
the  fluorescence,  thermally  assisted  fluorescence,  and  the  excita 
tion  pulse  allow,  evaluation  of  this  effect,  but  this  must  be 
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verified  for  significant  changes  in  flame  conditions;  and 

(ii)  if  the  ratio  of  fluorescence  (directly  pumped  by  laser  source) 

to  thermally-assisted  fluorescence  radiances  is  employed  with 

a  time  integrated  technique, then  further  assumptions  must  be 

made  in  order  to  determine  whether  E  >>  Es  ;  in  which  case, 

v0  v0 

the  thermally  assisted  fluorescence  population  rate  cannot 
possibly  follow  the  laser  pulse  and  the  two  signals  are  not 
correlated  in  a  single  way.  As  pointed  out  in  item  (i), 
thorough  studies  arc  needed  to  investigate  the  effect  of  the 

fiovy 

technique  (integrated  or  temporally  resolved  peak^is  affected 
by  the  temporal  shape  of  the  signals. 

FINAL  COMMENT 

All  five  methods  can  be  used  to  obtain  spatial  flame  temperatures, 
i.c.,  temperatures  corresponding  to^  1  mm3  of  flame  gases.  Methods  I, 

II,  III,  and  IV  can  be  used  for  temporal  flame  temperatures  as  long  as  the 
two  laser  beams  are  properly  matched  and  aligned  and  the  timing  between 
the  2  laser  pulses  is  properly  monitored;  temporal  temperatures  corresp¬ 
onding  to  ^  1  ys  should  be  obtainable  by  this  approach.  Method  V  is  ex¬ 
perimentally  implemented  most  easily  since  only  one  lase.r  pulse  is  needed 
which  simplifies  greatly  the  measurement  process.  Methods  II,  III,  and  IV 
strictly  require  saturation  resulting  in  problems  with  use  of  many  commer¬ 
cial  dye  laser  systems  especially  in  the  UV  when  doubling  of  the  dye  laser 
output  is  needed  with  the  concomitant  loss  in  spectral  irradiance  (also, 
unfortunately,  the  saturation  spectral  irradiance  aE*  depends  upon  v3xc 
and  so  much  more  source  spectral  irradiance  is  needed  in  the  UV  than  in 
the  visible  region  for  excitation).  Method  I  strictly  requires  linearity, 

between  B  and  E  ,  whereas  in  Method  V,  in  some  cases  E  is  limited  to  near 
v  v  o 

or  ' -"low  saturation.  In  Methods  III  and  IV,  certain  critical  assumptions  regarding 


88 


the  rate  constants  must  be  valid  to  obtain  accurate  flame  temperatures, 
and  fluorescence  measurements  must  be  made  in  the  portion  of  the  fluores¬ 
cence  pulses  where  the  3-level  steady  state  saturation  limitation  applies. 


In  Method  II  and  IV  ,  ,  no  assumptions  concerning  rate  constants  are 

pC  3  K 

needed,  but  the  experimental  system  must  allow  measurement  on  the  peak 
(2  level  steady  state)  of  the  fluorescence  pulse.  Methods  I,  II,  III,  and 
IV  all  require  careful  beam  matching,  alignment,  and  homogeneity,  whereas , 
Method  V  is  independent  of  these  factors.  Methods  II,  III,  and  IV  do 
not  require  spectrometer  calibration  as  long  as  the  same  fluorescence  line 
is  measured  for  the  ratio,  whereas.  Methods  I  and  V  require  calibration  of 
the  fluorescence  spectrometer  system  and  Method  I  requires  an  addition  cal¬ 
ibration  of  the  source  excitation  spectral  irradiance.  All  of  the  methods 
should  have  comparable  S/N  ratios  for  equivalent  metastable  energies; 
whereas  Method  I  will  have  the  lowest  average  signals,  the  increase  in 
measurement  time  should  produce  similar  S/N  ratios  to  methods  II,  III, 

IV,  and  I  (S/N  estimations  for  various  conditions  have  been  previously 
7 

estimated  )•  All  of  these  methods  are  based  on  pseudo-continuum  excitation 
conditions,  and  in  many  commercial  dye  lasets  detuning  or  other  line  broad¬ 
ening  techniques  must  be  initiated  to  assure  this  (when  E  <<  Ef  , 

'  v0  V0 
AXlaser  should  be  2*  °-2  ’  In  manY  cases,  even  though  the  total  laser 

energy  output  may  remain  about  the  same  after  line  broadening,  the  laser 
spectral  irradiance  may  be  severely  reduced. 

Experimental  flame  temperatures  obtained  by  these  methods  for 
hydrogen  based  flames  will  be  given  in  a  later  manuscript. 
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Figure  IB. 


Energy  Level  Diagram  of  3-Level  Tl-Like  Atom  With  Activation 
(Solid  Lines)  and  Deactivation  (Dashed  Lines)  Rate  Con¬ 
stants  Shown  (see  text  for  definition) . 
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V  , 


Figure  2B.  Hypothetical  Fluorescence  Temporal  Waveform  Produced  With 
a  Rectangular  Excitation  Pulse  (tgx  long).  Demonstrating 
the  2-Level  'Peak  and  the  3-Level  Steady  State  Plateau. 
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Hypothetical  Fluorescence  Temporal  Waveform  Produced  By 
Excitation  First  With  2-*3  to  Saturate  the  3  1  Fluorescence 
and  Then  While  2-»-3  Excitation  Occurs,  Excitation  Occurs' 
Simultaneously  at  2-*-3  and  l-*3  to  Saturate  the  3-*l 
Fluorescence.  During  the  Excitation.  Time,  tgJ.,  only  the 
2'*3  Excitation  Pulse  is  "On".  During  the  Excitation  Time 
Between  t  „  and  t  „  Both  2^3  and  l"*^  Excitation  Pulses 


I 


Figure  4B.  Energy  Level  Diagram  of  Multi-Level  Atom,  e.g.,  T1 ,  with 
Activation  (Solid  Lines)  and  Deactivation  (  Dashed  Lines) 
Rate  Constants  Shown  (see  text  for  definitions).  Level  2 
is  Shown  but  is  Not  Important  to  the  Processes  Shown. 
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ABSTRACT 


The  theoretical  basis  of  a  new  technique  called  thermally  assisted 
atomic  fluorescence  for  measurement  of  local  spa t ia 1  -  temporal  flame  temper¬ 
atures  is  given.  In  this  method,  the  ratio  of  fluorescence  signals  from  a 
radiatively  excited  level  and  a  higher  col 1  is ional ly  excited  level  or  from 
two  higher  col  1 i s iona 1 ly  excited  levels  of  an  inorganic  probe,  such  as 
In  or  T1 ,  is  related  to  the  flame  temperature.  The  conditions  necessary 
for  the  measured  flame  temperature  to  be  identical  to  the  translational 
(Boltzmann)  flame  temperature  are  given.  The  experimental  system  and  condi 
tions  for  making  flame  temperature  measurements,  which  consists  basically 
of  a  single  pulsed  dye  laser  and  a  gated  detector,  are  described.  The 
assumption  of  steady  state  conditions  is  discussed  as  well  as  several  ano¬ 
malies  in  terms  of  fluorescence  from  certain  energy  levels  of  T1 .  The 
accuracy  and  precision  of  the  measured  flame  temperatures  for  several 
methane-air  flames  arc  discussed. 
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INTRODUCTION' 


Laser  induced  fluorescence  (LIT)  is  one  of  the  most  promising  techni¬ 
ques  for  combustion  diagnostics,  as  it  is  spatially  ( and  temporally)  precise, 
very  sensitive  and  capable  of  detecting  trace  species.  Although  LIF  is 
mainly  employed  for  species  concentration  determination  [1,2],  it  seems  that 
temperature  measurements  can  also  be  performed. 

A  great  deal  of  work  has  been  done  to  assess  the  applicability  of  the 
two-line  atomic  fluorescence  method  [5,4,5]  to  flame  diagnosis,  while  other 
possible  techniques,  employing  LIF,  have  been  theoretically  developed 
[6,7]  . 

A  different  approach  consists  of  measuring  the  fluorescence  from  colli- 
sionally  excited  states  following  radiative  excitation.  If  the  collisional 
rates  are  fast  enough  to  allow  a  new  Boltzmann  distribution  during  the 
laser  pulse,  from  the  thermally  assisted  fluorescence,  the  local  tempera¬ 
ture  can  be  measured. 

A  partial  Boltzmann  equilibrium,  corresponding  to  the  flame  temperature, 
has  been  found,  for  Na  atoms,  in  a  H^-O^-Ar  flame  as  a  result  of  two-photon 
laser  excitation  [8].  Assuming  a  modified  two-level  model,  a  relatively 
low  deviation  from  the  Boltzmann  distribution  was  calculated  for  the  rota¬ 
tional  levels  of  the  Oli  molecule  [9],  and  rotational  and  vibrational  temp¬ 
eratures,  of  the  0(1  molecule,  as  a  result  of  collisional  energy  transfer, 
have  been  recently  measured  in  methane-air  flames,  using  laser  induced 
fluorescence  [10,11],  However,  the  lack  of  knowledge  of  the  magnitude 
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of  quenching  rates,  the  non- Bol tzmann Lan  rotational  population  distribu¬ 
tion  for  the  Oil  molecules,  and  the  dependence  of  vibrational  transfer  upon 
the  rotational  level  [11]  indicate  more  studies  have  to  be  carried  out, 
especially  for  molecules.  These  difficulties  can  be  partially  overcome 
by  seeding  the  flames  with  inorganic  atoms  and  measuring  well-characterized 
thermally  assisted  fluorescence  lines.  In  this  way,  spatially  resolved 
temperature  measurements  can  be  performed,  and  a  better  understanding  of 
the  collisional  energy  redistribution  can  be  determined. 

THEORETICAL  CONSIDERATIONS 

Two  Level  System 

It  is  generally  assumed  that  thermodynamic  equilibrium  exists  in  atmos¬ 
pheric  combustion  flames;  however,  this  is  not  "rigorously  true",  because 
of  radiative  losses  always  present  in  flames.  By  solving  the  rate  equa¬ 
tions  for  the  steady  state  case  for  a  two  level  system  (as  indicated  in 
Fig.  1)  in  the  absence  of  an  external  radiation  source,  the  ratio  (n,/nj) j. 


[CONTINUE  WITH  PAGE  4] 


is  given  by: 


(^)  =  (- 
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n2 

—x — )  tjrU 

A21  nl  B 

^71 


(i) 


where  (n2/n1)g  is  the  population  ratio  as  given  by  the  Boltzmann  relation¬ 
ship  and  A21  and  are  the  radiative  and  quenching  rates  from  level  2  to 
level  1. 

Hence,  from  the  ratio  (^/n^),  an  "optical”  temperature,  T ,  different 
from  the  "true"  translational  temperature,  is  given  by: 


£n(F 


'21 


21 


2 1 


(2) 


where  AE.^  is  the  energy  difference  between  the  levels  1  and  2,  and  k  is 
the  Boltzmann  constant.  This  well-known  deviation  is  normally  negligible 
because  <K  ^21’  ^ut  * t  *s  use^ul  to  remember  that  whatever  optical  dia¬ 

gnostic  technique  is  used  to  measure  the  ratio  n 2^n\>  tBe  resulting  temper¬ 
ature  will  be  the  "optical"  temperature,  corresponding  to  the  actual  pop¬ 
ulation  distribution  which  may  be  different  (frequently  very  little  differ¬ 
ent)  from  the  Boltzmann  distribution  obtained  by  the  translational  temperature. 

Three- level  system 

The  above  considerations  are  further  complicated  by  the  fact  that  for 
most  atomic  or  molecular  systems,  the  two- level  model  is  inadequate.  A 
three-level  system,  which  is  a  more  realistic  model  (sec  Fig.  2)  can  also  be 
solved  assuming  steady  state  conditions.  In  addition,  if  there  is  no 
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is  then 


radiative  excitation  by  an  external  source,  the  ratio  (n7/n^)jjj 
pi  veil  by: 


III 


k12^k32  *  k31  *  A31^  +  kl 3  k32 _ 

'(k2i  +  A21^k32  +  k31  +  A31^  +  k23^k31  +  A31^ 


(3) 


From  this  expression,  it  can  be  seen  the  (n2/n^)jjj  ratio  is  quite  different 
from  the  Boltzmann  ratio.  However,  if  it  has  been  previously  accepted  that 
the  n2/n.  ratio  is  Boltzmannian ,  then  one  is  forced  to  accept: 


n?  n 

^n^III  %  (-n~)II  ^ 


k12  =  £2 
S  8  "  *n  ”  8i 


ae2i 

exp  ( -pp=-i)  ; 
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(4) 


Therefore,  the  validity  of  the  Boltzmannian  equation  (4)  forces  the  assump¬ 
tion  that  =  k2„  =  0  and  A21  <<  k21,  or  that  the  radiative  losses  are 

negligible  (A21  <<  k21,  A31  <<  k31,  A32  <<  k32). 

Considering  now  the  case  in  which  an  external  source  of  spectral  irrad- 
iancc  E 

\ 

(refer  to  Boutilier,  et  al  [12])  by 


irradiates  a  three  level  system,  the  ratio  n-,/n7  is  given 
v12  ^  L 
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Expression  (5)  can  be  rearranged  as: 
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Also , 
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k^7  is  valid  in  most  three  level  cases  and 
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Then  it  can  be  readily  shown  that: 


where : 


n, 

o 


=  ( 


k23^ 

& 


B12  E12 
jqr  c 


2 1  - 


f  L*  *■  \ 

exP  C - ^ 


21 

*21 


B12  E12 

*77“ 


(i 


3 1  -j 

Y77> 

jl 


A3 1 

*77 


C  =  1 


‘23 


F7TY 


13  21 


( k !  3 


e12)  =  1 


k.. 

i 


;  3 7*1 2 


(k 


13 


‘12 


B12  E12^ 


101 


Several  limiting  cases  of  expression  8  are  now  possible.  If  no  external 
source  is  present,  the  ratio  n,/n2  is: 
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31 
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Once  again  the  two  levels  (3  and  2)  are  in  a  Boltzmann  equilibrium  only 

if  A2i  <<  k2^  and  A,^  <<  k.^,  i.e.,  the  radiative  losses  must  be  negligible. 

However  if  radiative  losses  are  negligible  (A's  <<  k's)  and  an  external  sourer 

of  spectral  irradiance  E  is  present,  then  the  ratio  n,/n7  becomes: 
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and  upon  rearrangement  and  substitution  for  C 
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Assuming  E^/c  >>  k^7  (a  condition  generally  achieved  with  laser 

i rradiances ) ,  the  ratio  is  given  by: 


7  , 


102 


(1 


n 

n 


3 

2 


■H 


(J^H 

k32 


exp  (- 


] : 


TT 

"T 


32 


32 


and  assuming  common  flame  temperatures  (>2000  K)  and  atomic  energy  levels 
(£l  eV) 


k23. 


=  (t7^)(t7 


'32 


- +~T — ^ 

31  K32 


(1 


Now  if  radiative  losses  are  not  neglibible,  i.e.,  A21  and  A,^  cannot 


be  disregarded  then  the  ratio  n,/n2  is 
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K12  K31  K31  K31  K12  K21 


®12  ^12^C 


"12, 


^32  ®12  E12/C 


'12 


A,.,  k,? 

f1  +  ir1  +  ir^  +  t1  +  +  f 

K31  K31  K31 


k32  ^  k32  k13  .  A31. 
31  ^  ^ 


Assuming  an  intense  external  source,  i.e.,  7  E^/c  >>  k12,  then: 


k7, 

2  j 


^21  k3? 

expC-fr)  + 


31 


F7T  A, 


'32  f'3 1  *  k31  +  k32 


'31 


d‘ 


(1 


and  assuming  common  flame  and  atom  conditions: 
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'32 


T7  '  '•p  +  X  +  p  ■> 

k32  31  a31  k32 


(16) 


Therefore  in  all  the  above  cases,  the  ratio  n./^  is  not  given  by  the  simple 
Boltzmann  relationship  of  and  only  approaches  it  in  a  limiting 

situation.  The  general  behavior  of  the  ratio  n^/n-2  qualitatively  shown 
in  Fig.  3. 

The  same  result  as  given  in  equation  (16)  can  be  more  easily  found  from 
equation  (5)  considering  =  0,  i.e.,  the  population  of  level  3  is  obtain¬ 

ed  only  by  collisional  coupling  with  level  2.  However,  the  assumption  k^,  = 

0  can  lead  to  some  misunderstanding.  Equation  (16)  is  independent  of  the 

E  value;  however,  for  low  irradiances,  the  assumption  k. ,  =  0  means  that 
v12 

n2k23  >:>  n1^13  and  ^32  >5>  ^31*  In  effect,  for  low  external  source  irrad¬ 
iances,  equation  (16)  must  reduce  to  the  Boltzmann  distribution. 

It  is  therefore  possible  to  determine  the  temperature  of  a  system  by 
measuring  the  ratio  of  two  fluorescence  lines,  one  coming  from  the  radiatively 

excited  level  (resonant  or  non-resonant,  Bc  ),  and  the  other  one  coming 

^  2 1  / 1 2 

from  an  upper  thermally  excited  level  (such  as  Br  ),  that  is: 

F  3 1 /  12 


BP  . 

*■31/12  _  ,v31,  ,*31wg3,  , 

B -  "  vtr  ’ >  (5 — J  (.5— )  (v - 

F  2 1  / 1 2  21  21  §2  31 


k32  ,  ,  AE32, 

+  A, ,  +’  exp  “FT"^ 


'31 


32 


(17) 


and  so 


T  = 


1E32/k 


*  tn(E 


'32 


21  8: 


31 


+  A 


^21  0  F  7  1  /  1 

*  ln(xrr)  * 

ol  0-  31  F31/12 


(18) 


104 


t  . 


This  technique  is  similar  to  the  classical  emission  two-line  ratio 
technique  for  temperature  measurement,  except  that  it  is  a  local  techni¬ 
que.  Of  course,  using  only  one  laser  excitation  line  and  observing  the 
two  fluorescence  lines  with  two  optical  detectors,  it  is  possible  to 
measure  spatially  and  temporally  resolved  temperatures. 

The  disadvantages  of  this  technique  are  the  need  to  have  accurate  A- 
values  (at  least  accurate  ratios)  and  to  have  an  accurate  calibration  to 
account  for  the  spectral  response  of  the  optical  set-up.  However,  this 
problem  is  minimized  by  using  a  calibrated  lamp.  The  value  of  the  ratio 
k,7/(k„7  +  A,,  +  k-.)  is  needed;  however,  for  common  flame  conditions  and 
with  seeding  atoms  like  Tl,  where  AEl2  a,3  eV  and  AE73  ^1  eV,  the  correction 
factor  in  equation  (16)  is  negligible,  i.e.,  is  essentially  unity. 


Four-level  system 

Now  consider  a  four  level  system,  where  the  4th  level  is  any  one  of 
the  levels  shown  as  a  manifold  of  levels  in  Figure  4.  Also,  let 
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The  rate  equations  can  then  be  written  as: 


d  n-, 

anr  =  VBi2  — 


Ev  (t) 
v12 


+  k12]  +  n3  A32  +  n4  A42  -  n2[A21  +  k23  +  k24] 


g-T-  "  n2  k2 5  +  nl  kl 3  +  n4  A4 3  '  n3[A32  +  A31  +  A34 


a- T"  “  n3  A34  +  n2  k24  +  nl  k14  '  n4 [ A4 1  +  A42  +  A43 


(20) 


nl  +  n2  +  n3  +  n4  =  nT 


With  the  assumption  k13  =  k14  =  0,  the  steady  state  solutions,  for  the 
ratios  n3/n2  and  n3/n4  are  given  by: 


n2 


k9_ 

(t4^)  (- 


32  1 


1  +  ( 


24, 

A? 


A4  3 


23 


+  /T* 

41  A42 


+  A x 
A4  3 


"*31 

aT 


A34 

+  jp-  a  - 

32  A32 


“A5 
A4  3 


A* 

A41 


(21) 


n_3 

n4 


‘23 


( A4 1  +  A4  2 ) 
,'43w  *17  ^23  +  k24) 


A 


7KWZ  +  aTTT 
,  .  -  ■  (  32  ol) 

*34  ^k23  +  k 


'24 


T 


(22) 


With  the  assumption  that  radiative  losses  arc  negligible,  i.e.,  A  <<  k, 
then 
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Now,  if  the  following  condition  is  va  lid  k31  <<  k32  +  k34 [k42^k42  +  k43)]> 
than  a  Boltzmann  distribution  occurs  between  levels  3  and  2. 

Remembering  that  level  4  can  be  of  any  level  in  a  manifold  of  levels 
(see  Fig.  4),  it  seems  that  the  general  conclusion  would  be  that,  apart 
from  the  radiative  losses  under  external  source  excitation,  in  a  steady 
state  case,  all  levels  above  the  excitation  level  are  always  in  Boltzmann 
equilibrium.  However  if  in  the  absence  of  external  source  excitation, 
levels  3  and  4  are  mainly  populated  via  the  coupling  k^2  +  ^3  anc* 
k12  +  k24  (assumPti°n  k^3  =  k ^  =  0) ,  then  a  Boltzmann  equilibrium  will  be 
also  between  these  levels  and  the  excited  level  2.  In  the  other  case  in 
which  the  main  thermal  population  of  levels  3  and  4  comes  from  level  1, 
through  k^3  and  k14>  laser  excitation  perturbs  the  distribution  and  the 
ratios  between  the  highly  excited  levels  and  the  radiatively  excited  level 
is  not  Boltzmannian .  However,  in  most  situations,  any  correction  made 
will  be  small. 

In  a  multi-level  system,  the  temperature  can  also  be  evaluated  using 

the  expression  (18),  with  or  without  the  correction  factor  that  can  be 

derived  from  Equation  (25),  but  now  there  is  the  possibility  of  measuring 

two  thermally  assisted  fluorescence  lines  (Bn  ,  Bc  ),  coming  from 

j/12  hi/12 

different  energy  levels  (Ej ,  Fh),  then  the  temperature  can  be  obtained  by 
the  expression: 


T 


V 


A .  g  .  X  .  F  •  /  i  -> 

£n  (-,1-  -)  +  SLn^)  +  £n(^ - ~) 


v  A .  g. 

J  h) 


l 


F; 


i/12 


(26) 


Another  possibility  is  to  consider  many  thermally  assisted  lines,  coming 
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from  different  energy  levels.  If  it  is  assumed  that  during  the  time  of  the 

laser  pulse,  the  system  reaches  a  new  thermodynamic  equilibrium,  from  a 

set  of  different  thermally  assisted  lines,  the  temperature  can  be  evaluated 

from  the  slope  of  a  plot  of  Sln(Bc  /  \> .  A.  g.)  versus  AE./k. 

hj/12  J  J  3  3 

EXPERIMENTAL 

The  apparatus  used  to  test  the  possibilities  of  the  thermally  assisted 
fluorescence  (THAF)  technique  is  reported  in  Fig.  5.  A  Nitrogen  laser 
(Molectron  UV-14)  pumped  dye  laser  (Molectron  DL-400)  is  used  for  this 
study.  The  optical  triggering  is  achieved  by  a  quartz  plate  and  a  photo¬ 
diode  (EGG-FND  100).  The  fluorescence  fluxes  are  collected  by  a  JY  H10 
monochromator  with  the  aid  of  a  quartz  lens  and  neutral  density  filters 
to  produce  approximately  the  same  fluorescence  signals  for  different  lines. 

A  Tektronix  1S1  (sampling  unit)  is  used  with  the  scan  facility  of  the  PAR- 
160  to  scan  a  fast  gate  (0.35  ns)  along  the  fluorescence  pulse.  The  sig¬ 
nals  are  visualized  by  a  Hewlett-Packard  1220  A  oscilloscope  and  a  MFE 
plotmatic  x-y  715  recorder. 

For  this  study,  there  was  the  necessity  of  a  detector  with  a  response 
time  fast  enough  to  follow  the  fluorescence  pulses.  This  was  achieved 
using  a  Hamamatsu  R928  squirrel -cage  photomultiplier,  supplied  with  high 
voltage  (-1400  V)  and  a  laboratory  constructed  voltage  divider  following 
the  suggestion  found  in  the  literature  for  nanosecond  studies  [13,14,15,16). 

With  the  apparatus  used,  for  transient  pulse  measurements,  rise  times 
less  than  1  ns  are  achieved,  but  unfortunately  ringing  appeared.  However, 
we  were  concerned  with  reaching  steady  state  conditions  and  therefore  with 
mainly  the  "rising  portion"  and  not  with  the  decay  portion  of  the  fluorescent 
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transient,  and  the  presence  of  a  ringing  could  be  ignored  when 
signals  at  the  peak  of  the  curve  where  steady  state  conditions  could  be 
assessed. 

One  of  the  aims  of  this  work  is  to  investigate  the  possibility  of 
measuring  the  fluorescence  from  highly  excited  levels  even  for  quenching 
flames,  and  so  a  methane-air  flame  is  chosen  for  the  experiments.  The 
gases  are  regulated  by  precision  rotameters  (Matheson)  and  the  metal  solu¬ 
tions  are  seeded  into  the  flame  by  pneumatic  nebulization  with  a  concentric 
glass  nebulizer  (T-230-A3  J.E.  Meinhard) .  1000  ppm  solution  of  In  and 
T1  are  used  to  investigate  four  flames  with  composition  given  in  Table  I. 


RESULTS  AND  DISCUSSION 

For  these  investigations,  steady  state  conditions  must  be  reached 

during  the  laser  pulse.  A  preliminary  investigation  on  T1  in  a  methane- 

air  flame  indicated  that  all  thermally  assisted  fluorescence  lines  with 

upper  levels  within  1.5  eV  above  the  radiatively  excited  level  have  their 

fluorescence  peak  at  the  same  time,  roughly  1  ns  after  the  fluorescence 

(Fig.  6) 

peak  from  the  radiatively  excited  leveM  The  redistribution  of  the  popula¬ 
tion  among  the  upper  levels  seems  to  happen  approximately  at  the  same 
time,  and  steady  state  conditions  seem  to  be  reached  just  before  the  end 
of  the  laser  pulse.  Moreover,  the  good  agreement  found  in  the  temperature 
measurements  confirmed  that  the  laser  pulses  used  are  sufficiently  long 
to  reach  steady  state  conditions  at  least  in  the  methane-air  flame. 

The  first  element  used  as  a  thcrmometric  species  was  In;  in  Fig.  7 
a  simplified  energy  scheme  with  the  transitions  investigated  is  given. 

The  transition  between  the  6s  level  and  the  ground  state  is  pumped  by  the 

o 

dye  laser  and  the  fluorescence  of  the  resonant  line  (4101.76  A),  the 
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3258-3256  A  doublet  from  the  5d  levels  and  the  3039  A  line  are  measured. 


The  other  possible  lines  of  In  happened  to  interfere  with  each  other  and 
are  not  measured.  In  Tig.  8  a  typical  fluorescence  pulse  from  the  5d 
level  in  given  as  recorded  with  our  apparatus. 

The  temporal  shape  of  the  signal  is  well-characterized  with  a  distinct 
rise  time  and  a  well-defined  decay  curve.  The  temperature  found  in  the 
middle  of  two  slightly  different  flames  are  reported  in  Table  II.  The 
good  agreement  between  the  temperature  measured  suggests  that  the  preci¬ 
sion  of  the  technique  is  good  and  the  errors  can  be  less  than  a  few 
percent . 

More  measurements  were  performed  using  T1  as  the  seeding  element.  The 
energy  levels  and  the  fluorescence  lines  observed  for  the  first  set  of 
experiments  with  T1  are  reported  in  Fig.  9.  The  7s  level  from  the  ground 

o 

state  is  excited  and  fluorescence  at  the  strong  resonant  (3776  A)  and 
stokes  line  (5350  A),  and  ttie  3519-3529  A  doublet,  from  the  6d  thermally 
excited  levels  are  all  observed. 

The  results  and  the  temperatures  are  in  Table  II.  One  consideration 
is  evident  at  this  point.  The  temperatures  calculated  using  the  intensities 
of  the  3776  A  line  are  about  100  K  less  than  those  using  the  5350  A  line. 
Comments  on  this  will  be  given  later  in  the  paper.  One  of  the  attributes 
of  this  technique  is  that  it  is  possible  to  pump  highly  excited  levels  where 
several  energy  levels  within  £1.5  eV  arc  present. 

In  a  second  set  of  measurements  on  T1  the  6d  levels  are  pumped  by 
frequency  doubling  the  dye  laser  output;  fluorescence  lines  coming  from 
the  6d ,  8s  and  7d  levels,  could  all  be  observed  (see  Fig.  10).  The 
results  are  summarized  in  Table  II.  Considerable  discrepancies  arc  found 

o  o 

for  temperature  cn 1  oil  a t i ons when  considering  the  2767  A  or  the  3525  A 


1 1 1 


t'  1  uo  re  scene  o  linos  from  the  same  (id  level. 

o 

Because  the  377o,  S.i50  ,  2767  and  3S2.ri  A  fluorescence  lines  give  the 
same  results  (temperatures  1  with  either  the  boxcar  averager  with  a  15  ns 
gate  width  or  with  the  Tektronix  1S1  sampling  unit,  appreciable  ionization 
is  not  present  when  using  a  5  ns  laser  pulse  [17]. 

With  integrated  measurement,  i.e.,  the  boxcar  averager  approach,  ex¬ 


citing  the  7s  level  with  the  3776  A  line  fluorescence  lines  coming  from 
the  7s,  6d,  8s  and  7d  levels  are  observed.  The  latter  is  1.9  eV  above 


Bf-X 


the  excitation  level.  Plots  of  Hn (— )  versus  level  energy  are  reported 
in  Fig.  10  for  both  the  temporally  resolved  (1S1  plug-in)  and  integrated 
(boxcar  averager)  techniques.  It  is  apparent  that  the  fluorescence  lines 


at  5350,  3525  and  3229  A  result  in  a  straight  line,  corresponding  to  a 
1626  K  temperature  and  showing  a  Boltzmannian  distribution  for  the  popula- 

o 

tion.  On  the  other  hand;  the  fluorescence  lines  at  3776,  2767,  and  2920  A 
deviate  from  the  straight  line  portion.  The  optical  system  was  recalibrated, 
but  no  particular  differences  were  found.  The  only  explanation  at  this 
point  is  that  the  A-vaules  reported  in  the  literature  [18]  are  not  correct. 


CONIC  1;US  IONS 

We  have  performed  a  theoretical  analysis  of  the  population  ratio  for 
thermally  excited  levels  in  a  multilevel  scheme  and  steady  state  conditions 
and  we  have  shown  the  possibility  of  using  thermally  assisted  lines  to 
measure  the  temperature  of  the  flames.  We  have  also  performed  an  experi¬ 
mental  study  on  the  applicability  of  this  technique  on  a  methane-air  flame 
and  we  have  shown  that  using  particular  lines  accurate  temperature  measure¬ 
ment  can  be  performed.  However,  more  studies  are  needed  to  define  the 
correct  A'values  for  the  atomic  transitions  and  to  show  the  applicability 
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of  this  technique  in  different  experimental  condi 
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TABLE  I 


Gas  Composition  For  the  Four  Flame  Investigated 


Flame  Methane  Air 

_ (E/m) _ (E/m) 


1 

.50 

6.1 

2 

.51 

6.1 

3 

.50 

7.1 

4 

.50 

6.6 

1 
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TABLE  I  I  . _ S 1 1 inmary  of  the  Temperatures  Measured 

in  Throe  Different  Flames 


Flame 

Dine  Observed 

Level  Energy 

Temperature 

(A) 

(cm"1) 

00 

4102. 

24  373 

1725-1699 

3257  . 

32900 

4102 

24  37  3 

1785 

3257 

32900 

4102 

24  373 

1769 

3039 

32892 

3776 

26478 

1465-1426 

3525 

36159 

5350 

26478 

1565-1539 

3525 

36159 

2768 

36118 

743 

3229 

38746 

3525 

36159 

1442 

3229 

38746 

3229 

38746 

2408 

2920 

42030 

3525 

36159 

1843 

2920 

42050 

2768 

36118 

1199 

2920 

42030 

a : 

I  n 

as 

thermometr ic 

b: 

T 1 

as 

thermometric 

c : 

T1 

as 

thermometric 

species.  Excitation 
species.  Excitation 
species.  Excitation 


1 ine  at  4101.76  A 
line  at  3775.7  A 
line  at  2767.8  A 
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Figure  2C.  Three-Level  System  Under  Spectral  Irradiances  E 

Collisional  Excitation  of  the  3rd  Level.  12 

Radiative  and  Collisional  Rate  Constants  Are  Shown 
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and 


m 


Figure  4C.  Four-Level  System  Under  Spectral  Irradiance  Ey  and  with 
Radiative  and  Collisional  Coupling  Among  12 

the  Levels.  Significant  Ratios  of  the  Population  are  Given 
in  the  Text. 
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Figure  5C.  Experimental  Set-Up  for  the  Thermally  Assisted  Fluorescence  Measurements- 
N2L  -  Nitrogen  Laser;  DL  -  Dye  Laser;  PD  -  Photodiode;  B  -  Boxcar;  . . 

rF„  “  Flame;  ND  -  Neutral  Density;  M  -  Monochromator;  PM  -  Photomultiplier; 


in 

u 

o 


•H 

c 

3 

00 

C 

•H 


cu 

E 

m 

X 

•H 


a 

o 

u 

•M 

X 

0) 

H 


>s 

r-l 

P. 

ex  . 
p  X 
n  o) 
T3 
S-  u 
<D  O 

£  u 

O  ID 

x  a: 

<1)  I 

oo 

cO  OS 
M 

H  •  - 

O  0) 

>  n. 

o 
c  u 

00  to 

H  o 

n  i— ) 


u 
>  to 

n  o 


Figure  6C.  Temporal  Resolved  Fluorescence 0Pulses  of  T1 .  The 

Thermally  Assisted  Line  (3525  A)  has  its  Peak  Roughly 
1  ns  After  the  Fluorescence  Maxim  from  the  Radiatively 
Excited  Level  (3776  A) . 


Scheme  of  T1  Showing  the  Transitions  Observed  (arrows).  The  Energy 
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ABSTRACT 


The  broadening  of  the  absorption  profile  in  hydrogen-based  flames  diluted 
with  argon  and  nitrogen  is  discussed  theoretically  and  demonstrated  experimentally 
with  a  pulsed  tunable  dye  laser.  This  broadening  is  evaluated  from  the  half¬ 
width  of  the  profile  obtained  by  scanning  the  laser  beam  through  the  atoms  in 
the  flame  while  monitoring  the  resulting  fluorescence  with  a  high  luminosity 
monochromator,  i.e.,  by  observing  a  fluorescence  excitation  profile.  Results 
are  given  for  the  elements  Ca,  Sr,  Na  and  In.  It  is  shown  that  the  halfwidth 
of  the  atomic  profile  depends  approximately  upon  the  square  root  of  the  log  of 
the  laser  irradiance.  This  dependence  stems  from  a  theoretical  treatment  based 
upon  the  interaction  of  a  gaussian  laser  profile  and  a  gaussian  atomic  profile. 

It  is  also  shown  that  if  the  laser  spectral  bandwidth  exceeds  by ~  5-10  times 
the  absorption  halfwidth,  the  fluorescence  excitation  profile  provides  a  simple 
means  of  a  reliable  evaluation  of  the  laser  bandwidth. 
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INTRODUCTION 


The  high  irradiance  provided  by  the  pulsed  tunable  dye  laser  is  capable 

of  saturating  both  single-photon  and  2-photon  transitions  of  atoms  in  flames 

■1-7 

at  atmospheric  pressure.  Among  other  effects,  such  as  the  attainement 

of  the  maximum  fluorescence  signal  and  its  relative  independence  upon  the 
quantum  efficiency  of  the  transition,  the  strong  irradiation  field  is  known 
to  be  responsible  for  a  broadening  of  the  absorption  line  profile.  This 

g 

broadening  is  called  saturation  broadening.  Hosch  and  Piepmeier,  Van 
4  9 

Dijk,  as  well  as  others  reported  an  experimental  verification  of 
saturation  broadening  by  performing  both  absorption  and  fluorescence  measure¬ 
ments  . 

The  aim  of  this  work  is  to  present  a  theoretical  treatment  of  the  saturation 
broadening  based  upon  the  rate  equations  approach  and  for  the  two  cases  of  a 
very  narrow  excitation  line  laser  coupled  with  a  Lorentzian  absorption  profile 
and  a  gaussian  laser  profile  coupled  with  a  gaussian  atom  profile.  The  theoretical 
predictions  will  be  compared  with  the  experimental  data  obtained  by  scanning 
the  laser  beam  throughout  the  atom  profile  in  different  flames  while  monitoring 
the  fluorescence  emitted  at  right  angles,  i.e.,  by  obtaining  a  fluorescence 

c-n  10,11 

excitation  profile. 

We  realize  that  our  experimental  conditions  are  not  applicable  to  the  first 
case  (line  source  and  Lorentzian  atom  profile).  Furthermore,  the  second  case 
(gaussian  laser  profile  and  gaussian  atom  profile)  is  certainly  an  approximation 
to  the  real  interaction  process.  Nevertheless,  we  felt  that  it  was  useful  to 
present  our  experimental  data  since  the  saturation  broadening  was  clearly  ob¬ 
served  in  all  the  cases  investigated. 
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THEORY 


The  treatment  given  here  follows  the  basic  discussion  on  saturation 
which  can  be  found  in  several  books,  chapters  and  articles.1  7  Generally, 
atoms  are  considered  to  be  dispersed  as  trace  constituents  in  a  gas  at  at¬ 
mospheric  pressure  and  characterized  by  a  2-level  system.  Coherence  effects 

(cooperative  phenomena)  are  neglected  because  dephasing,  coherence-interrupting 

12  13 

collisions  are  considered  to  be  fast  in  our  atmospheric  pressure  flames.  ' 


Case  A.  Line  source  and  Lorentzian  atom  profile 

In  this  case,  the  laser  is  assumed  to  be  a  very  narrow,  monochromatic  line 

source  and  the  atomic  absorption  profile  is  considered  to  be  homogeneously 

14  is 

broadened.  The  first  assumption  can  be  met  by  tunable  dye  lasers.  Also, 

the  laser  beam  is  assumed  to  be  spatially  uniform.  The  case  in  which  a  gaussian 
spatial  profile  is  assumed  for  the  radial  dependence  of  the  laser  intensity16,17 
is  discussed  in  the  Appendix. 

The  atom  profile  can  be  considered  homogeneous  if  the  following  conditions 

hold:  (i)  Doppler  broadening  is  negligible  compared  to  collisional  broadening; 

(ii)  in  the  case  of  combined  Doppler  and  collisional  broadening,  velocity  changing 

collisions  are  so  fast  that  atoms  cannot  be  considered  to  belong  to  any  particular 

Doppler  interval  during  the  time  of  interaction  with  the  laser  beam. 

Let  the  atoms  be  characterized  by  levels  1  and  2,  with  energy  difference  hVQ(J), 

statistical  weights  (dimensionless)  g.  and  g  and  number  densities  (m  3)  n1  and  n  . 

12  2 

The  total  atomic  population  in  this  case  is  nT=  +  n2"  In  tIle  steady-state 
limit,  we  can  obtain  from  a  simple  rate  equations  approach,7  the  following 
relationship: 
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(1) 


B12  (»<V/8X.ff1  9  <VV 


'“'(■“feK'v 


V  +  *21  *  *21 


where  ■-  .... 

=  Einstein  coefficient  of  induced  absorption,  J  Xm^s  '''in;  =  Einstein 

coefficient  of  spontaneous  emission,  s  *  k^  =  collisional  (quenching)  rate 


constant 


,  s  X  p(X^)  -  I  p^(X)dX=p^  dX^  ,  integral  energy  density  of  the 

J  SX, 


£.  -3 

laser,  integrated  over  the  effective  width  of  the  laser  line,  6X^  ,  in  Jm  ; 

gJ<X  -X  )  =  absorption  shape  function,  given  by  the  Lorentzian  dispersion  formula 

6X_ 
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6X  =  full  width  at  half  maximum  of  the  shape  function,  m  ; 
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6X  ^  =  effective  width  of  the  absorption  profile,  m  ; 
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dimensionless. 


Equation  (l^  can  be  rewritten  in  terms  of  the  effective  lifetime,  T  (flA^H-  ^21'  ^ 
of  the  excited  level  as 


n2  B12{p(V/6Xeff}  ^VV  T 


(1  +  ~  )  B12{p(Xn)/<5X_cc}  ga(X0-Xj  r  +  1 
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If  the  atomic  density  is  low,  then  n can  be  expressed  in  terms  of  the 

-1  -2  -1. 

fluorescence  radiance,  B  (J  s  m  sr  ) ,  as 

F 


BF  n2  (~)  A21  hV 
4tt  O 


(2) 


where  Z  (m)  is  the  homogeneous  depth  of  the  fluorescence  volume  in  the 
direction  of  observation. 


We  now  define  the  saturation  energy  density 
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pS<V 

eff 


A21  +  k21 


d+Ii)  B12 
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(  1  +  li)  B12T 
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d  +  ^)B12Y21 


(3) 


where  Y  =  {a  /(A  +  k  )}  is  the  quantum  efficiency  of  the  transition. 

4  X  «  X  4  1  4  X 

By  rearranging  Eq.^l^and  making  use  of  Eqs.  (2)  and  (3),  we  obtain  the  atomic 
fluorescence  radiance  the  following  expression: 


where  C  =  hV  (Z/4tt) 
0 


A 


21’ 


This  equation  shows  the  theoretical  dependence  of  the  wavelength  integrated 
fluoresence  radiance  upon  the  laser  energy  density  and  the  detuning  between  the 
peak  of  the  laser  profile,  acting  here  as  a  delta  function,  and  the  peak  of  the 
Lorentzian  atom  absorption  profile.  It  can  be  easily  shown  that  this  equation, 
in  the  low  intensity  limit,  reduces  to  the  well  known  linear  dependence  of  the 
fluorescence  radiance  upon  p(X^)  and  Y^  while,  in  the  high-intensity  limit,  n^ 
becomes  half  of  the  total  atomic  population  if  the  statistical  weights  are 
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the  same.  This  particular  outcome  holds,  of  course,  irrespective  of  the  line 


shape . 

The  ratio  between  braces  in  Eq.  ^  is  a  modified  Lorentzian  function,  which 
takes  into  account  the  effect  of  the  laser  field  upon  the  spectral  profile. 

The  low  intensity  Lorentzian  function  can  be  written  as 


6X. 


(5) 


,6XL.2  .  ,  .2 

(— )  +  a-\0) 


while  the  modified  Lorentzian  function  in  Eq.  ^  can  be  written  as  a  function  of 


the  laser  density  p  and  the  detuning  6  =  (X  -X  )  as  follows  : 

J6  0 


In  both  equations,  <5X  is  full  width  at  half  maximum  (FWHM)  of  the  unperturbed 

L 

Lorentzian  function.  Inspection  of  Eqs.  (5)  and  (6)  shows  that  the  FWHM  of 

the  modified  function,  here  called  6X  . ,  ^ .  ,  is  given  by 

excitation 


6X  ,2  6X  ,  2  6X  .  2  |P(V 

(  exc)  =  (  L)  +  {  L) 


pS(V 


or 


<5A 


exc 


P(V 

pS(X£) 


(7a) 


(7b) 


This  equation  shows  that  the  width  of  the  absorption  profile  increases  with  the 
square  root  of  the  laser  density. 


(i)  Since  pS(A^)  depends  upon  the  quantum  efficiency  of  the  transition,  the 
broadening  effect  at  a  fixed  laser  density  will  be  greater  in  high  quantum 
efficiency  flames  as  compared  to  that  observed  in  low  quantum  efficiency  flames. 

(ii)  As  indicated  by  Eq.  ^7a^,  a  plot  of  the  experimentally  observed  (5X  c)^  vs. 

the  laser  density  measured  at  the  flame  should  result  in  a  straight  line.  From 
the  slope  of  this  line,  one  can  calculate  the  saturation  density  for  that  parti¬ 
cular  transition  while  from  the  intercept  one  can  calculate  the  low  intensity 
Lorentzian  linewidth. 

(iii)  At  the  limit  of  p(A.)  =  0,  6A  =  6A  ,  as  it  should  be. 

(iv)  Since  the  second  term  in  the  square  root  involves  a  ratio  between  the 
densities,  one  needs  only  to  measure  the  laser  power  $,  at  the  flame. 

Equation (V^can  also  be  derived  from  Eqs.  (4)  and  (6)  in  the  following  manner. 
When  the  detuning  in  zero,  the  fluorescence  signal  is  peaked  at  a  value  given  by 
the  following  expression: 


when  the  detuning  is  very  large,  the  fluorescence  signal  approaches  zero. 

+ 

Therefore,  we  can  find  the  detuning,  6,  ,  at  which  the  peak  value  of  the  function 

is  halved,  for  both  positive  and  negative  detunings.  We  then  have 
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which  gives 


<5A 

exc 


+  (6)p 


p(A  ) 

1  + - — - 

0s ‘V 


and  is  seen  to  be  identical  with  Eq.  (7) . 
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Case  b.  Gaussian  laser  profile  and  gaussian  atom  profile 

This  case  can  be  considered  to  be  closer  to  our  experimental  set-up  when 
the  pulsed  laser  beam  is  not  very  narrow  and  the  atom  profile  in  the  flame  is 
described  by  a  Voigt  profile.  Here,  the  laser  is  assumed  to  be  gaussian  so 
that  its  spectral  energy  density  (J  m  3  m  S  is  given  by  the  following  expression 


Px (A) =  p 


2/  In 


TT  6A, 


exp 


{- 


2/&n  2 

6A„ 


(A-A^r 


(11) 


where  p  is  the  integrated  energy  density  (J  m  "S  and  6A^(m)  is  the  FWHM  of  the 
laser  spectral  profile.  Accordingly,  the  atom  profile  is  given  by  the  following 
dimensional  (m  S  shape  function: 


a  ,  ,  2  /  In  2  r 

gA  *  ~  o  =  — — ~ -  exp 


2/  Sin 


/  ir 


5A 


6A 


(X-X0))‘ 


(12) 


where  <5A  is  now  the  FWHM  of  the  atom  spectral  profile.  As  stated  before,  the 
atoms  here  are  not  considered  to  be  grouped  in  Doppler  intervals  (which  would 
cause  hole  burning  for  a  spectrally  narrow,  saturating  laser  beam)  because  of 
the  very  effective  cross-relaxation  taking  place  in  the  flame.  The  interaction 
of  the  laser  and  the  atoms  in  Eq.  (-’)  i  s  now  given  by  the  convolution  integral  of 
both  profiles.  Since  the  convolution  of  two  gaussian  functions  is  still  a 
gaussian  function,  we  obtain 


t 
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(13) 
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s  s  -1-2 

where  4>  and  <J>  correspond  to  p  and  p  but  are  fluxes,  in  J  s  m 

2 

From  Ect.  (18),  one  can  see  that,  by  plotting  vs.  £n$,  a  straight 

» 

line  is  obtained.  The  slope  of  this  line  gives  6A  ,  while  the  intercept  gives 
.  s 

. 

EXPERIMENTAL 

IS 

The  experimental  set-up  used  is  described  in  detail  elsewhere.  The 

tunable  dye  laser  (DL-400,  Molectron,  Sunnivale,  CA.)  pumped  by  a  nitrogen  laser 
(UV-14,  Molectron,  Sunnivale,  CA)  is  directed  into  the  flame  by  means  of  two 
plane  mirrors  after  passing  through  two  spatial  filters  consisting  of  shielded 
iris  diaphragms  (Edmund  Scientific,  Barrington,  N.  J.).  These  filters  passed 
only  the  central  portion  of  the  beam  so  to  improve  its  uniformity.  No  lenses 
were  placed  between  the  iris  and  the  flame.  The  flame  used  was  an  oxygen-hydrogen 
mixture  diluted  with  argon  or  nitrogen,  and  supported  by  a  capillary  burner. 

For  these  measurements, the  flame  was  surrounded  by  a  gas  sheath  (Ar  or  N^) 
but  was  not  protected  by  a  similar  analyte-free  flame  burning  at  the  same  composi¬ 
tion  and  temperature.  All  measurements  were  taken  at  approximately  1  cm  above 
the  primary  reaction  zone.  The  standard  solutions  for  all  elements  (Ca,  Sr, 

Na  and  In)  were  made  from  reagent  grade  chemicals.  The  concentration  was  chosen 
to  be  very  low  (in  most  cases  1  yg/mL)  to  avoid  any  self-absorption,  pre-filter 
or  post-filter  effects. 

The  fluorescence  was  collected  at  right  angle  by  an  optical  system  consisting 
of  two  spherical  quartz  lenses  and  a  rectangular  aperture.  The  flame  was  imaged 
at  unit  magnification  onto  this  aperture  and  after  that,  with  the  same  magnifica¬ 
tion,  onto  the  entrance  slit  of  a  high  luminosity  monochromator  (Jobin  Yvon, 
model  H-10,  10-cm  focal  length,  f-3.6)  whose  slit  width  was  set  for  all  elements 
at  0.5  mm.  The  signal  from  a  photomultiplier  wired  for  fast  pulse  high  current 
work  arid  operated  at  -1000  V  was  processed  by  a  dual  channel  boxcar  integrator 
(model  102-164,  PAR,  Princeton,  N.  J.).  In  order  to  sample  only  a  small  portion 


1  38 


of  the  fluorescence-  waveform,  which  in  our  case  is  the  result  of  the  saturation 
process  and  by  the  photomultiplier  temporal  characteristic's,^  ,  we  have  used 
a  75-ps  risetine  sampling  head  ( PAR  Model  163  sampled  integrator).  The  mono¬ 
chromator  was  centered  at  the  peak  of  the  fluorescence  line  and  the  laser  out¬ 
put  was  slowly  scanned  across  the  absorption  profile.  These  measurements  were 
then  repeated  by  decreasing  the  laser  irradiance  with  calibrated  neutral  density 
filters.  Particular  care  was  taken  to  minimize  spurious  reflections  and  scatter 
of  the  laser  into  the  monochromator. 

The  peak  power  of  Lhe  laser  at  the  particular  transition  used  was  measured 
with  a  calibrated  photodiode  (Model  F4000,  ITT,  Fort  Wayne,  Indiana)  coupled 
directly  to  a  scope  (Type  454,  Tektronix,  Portland,  Oregon). 

The  laser  spectral  bandwidth  was  measured  directly  with  a  good  resolution 
grating  monochromator  (Jobin  Yvon,  Model  HR-1000,  1-m  focal  length,  f-5.4  ). 

The  spectral  slit  function  of  such  monochromator,  evaluated  by  scanning  a  very 

O 

low  pressure  mercury'  pen  light  in  front  of  it,  was  found  to  be  0.12  A. 

RESULTS  AND  DISCUSSION 

Figures  1-4  show  clearly  the  broadening  effect  on  the  line  profile  as  the 
laser  irradiance  on  the  atoms  increases.  The  broadening  in  the  halfwidth  ranges 
from  a  factor  of  1.3  for  In  to  5.4  for  Sr.  Obviously,  as  expected  and  stated 
before,  this  broadening  cannot  be  accounted  for  by  Ea.  (7b)  which  was  derived 
with  the  assumption  of  a  monochromatic  laser  and  a  Lorentzian  atom  profile. 
Indeed,  since  in  some  cases  (Sr  and  Ca),  the  laser  power  exceeded  the  saturation 
power  by  approximately  a  factor  of  100,  Eq.  (7b)  would  have  predicted  a  10- 
fold  broadening  of  the  profile.  Except  for  the  case  of  indium  (Fig.  1)  for  which 
a  15  ns  boxcar  gate  width  was  used  to  average  the  signal,  the  profiles  have  been 
obtained  with  the  75  ps  risetime  sampling  head.  This  explains  the  different 


noise  levels  shown  in  the  figures. 


Irrespective  of  the  theoretical  assumptions  made,  both  Eq .  (7)  and, 

(lg\  predict  a  larger  broadening  effect  with  high  quantum  efficiency  flames  as 

compared  with  that  observed  with  low  quantum  efficiency  flames.  Despite  the 

fact  that  our  flames  are  not  shielded  (apart  from  an  outer  gas  sheath) ,  the 

replacement  of  argon  as  diluent  by  nitrogen  should  change  considerably  the 

quantum  efficiency.  Thus,  different  broadenings  are  expected  for  the  same 

laser  power  but  for  the  two  flames.  The  experimental  results,  collected  in 

Table  I,  do  indeed  show  the  correct  trend,  with  the  exception  of  strontium. 

These  values  are  therefore  indicative  of  the  different  quantum  efficiency 

and,  in  turn,  of  the  different  saturation  powers  in  the  two  flames.  This 

19 

matter  is  discussed  in  greater  details  m  another  paper. 

In  an  attempt  to  compare  the  experimental  data  with  the  theoretical  pre¬ 
dictions  given  by  Eg.  (16) ,  we  have  calculated  for  both  flames  the  ratio 

between  the  maximum  value  of  6A  (corresponding  to  that  obtained  at  full 

exc 

laser  power)  and  the  minimum  value  of  <5^exc  (corresponding  to  that  obtained 

at  p  «  p"  ,  i.e.,  when  6A  c  ~  "  The  saturation  power  was  evaluated 

19 

for  each  element  from  the  experimental  saturation  curve.  in  Table  II,  we  compare 

this  ratio  with  those  obtained  experimentally,  where  the  values  for  ( <5 A  ) 

exc  mm 

are  those  obtained  at  very  low  laser  powers  (see  also  Table  III) .  Apart  from 
some  unexplained  discrepancies  outside  the  range  of  the  experimental  errors 
(such  as  the  experimental  value  for  strontium  in  the  nitrogen-diluted  flame) , 
the  agreement  between  the  theory  and  the  experiment  can  be  considered  fairly 
satisfactory,  when  all  the  assumptions  made  at  the  beginning  are  properly  and 
critically  considered. 

In  the  flames  used  in  this  work,  the  halfwidth  of  the  atom  profile  is 

O 

expected  to  be  of  the  order  of  0.05  A  for  all  the  elements  considered.  Therefore, 
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fcii  ii  - 

the  convolution  obtained  by  scanning  the  laser  through  the  absorption  profile 
at  low  laser  powers  can  be  practically  entirely  attributed  to  the  laser  spectral 
halfwidth.  For  our  theory  to  be  consistent,  the  values  so  obtained  (fluorescence 
excitation  profile)  should  be  the  same  as  those  given  by  the  slope  of  the  plot 
of  Eq.  (16),  as  explained  before.  Table  III  shows  the  different  values 
obtained  with  these  two  methods.  In  addition,  the  values  obtained  by  measuring 
directly  the  laser  profile  with  a  medium  resolution  monochromator  are  also  shown. 
Again,  apart  from  some  differences,  the  average  value  for  the  two  flames  obtained 
with  both  techniques  (fluorescence  excitation  profile  and  saturation  broadening) 
closely  agrees  with  that  measured  directly.  We  therefore  feel  that,  whenever  the 
laser  bandwidth  exceeds  by  approximately  5  to  10  times  the  atom  profile,  the 
halfwidth  obtained  from  the  low  intensity  fluorescence  excitation  profile  provides 
a  reliable  estimate  of  the  laser  spectral  bandwidth.  Obviously,  this  method  fails 
when  the  laser  bandwidth  is  much  narrower  than  the  atom  profile  (since  in  this 
case  the  convolution  will  be  essentially  given  by  the  atom  profile)  and  loses 
its  significance  when  the  laser  bandwidth  is  much  larger  (>  20  times)  than  the 
atom  profile.  In  this  last  case,  it  would  be  indeed  much  simpler  to  scan  the 
laser  directly  through  a  small-medium  resolution  monochromator. 

CONCLUSIONS 

As  stated  at  the  beginning,  it  was  not  the  purpose  of  this  work  to  find  a 
theoretical  treatment  of  the  saturation  broadening  that  would  perfectly  fit  our 
experimental  conditions.  In  fact,  in  order  to  do  this,  the  exact  spectral  shape 
of  the  laser  pulse  and  the  atom  profile  have  to  be  known.  Nevertheless,  we 
feel  that  the  results  obtained  are  useful  and  can,  at  least  qualitatively,  be 


explained  by'  our  theoretical  approach. 


The  main  results  of  this  work  may  be  summarized  as  follows: 


(i)  As  has  been  known  theoretically  and  shown  experimentally  before, 
the  atomic  profile  broadens  when  the  laser  irradiance  is  such  that  saturation 
can  be  approached. 

(ii)  The  broadening  effect  is  larger  for  higher  quantum  efficiency  flames. 

(iii)  For  a  pulsed,  tunable  dye  laser  (without  spectral  narrowing  elements 
in  the  cavity) ,  pumped  by  a  nitrogen  laser  and  a  flame  at  atmospheric  pressure, 
the  broadening  depends  approximately  upon  the  square  root  of  the  log  of  the  laser 
irradiance . 

(iv)  The  halfwidth  of  the  fluorescence  excitation  profile  obtained  at  low 
laser  powers  can  indeed  be  taken  as  a  measurement  of  the  laser  spectral  band¬ 
width  if  this  bandwidth  is  approximately  5-10  times  larger  than  the  atom  profile. 
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APPENDIX 


Derivation  of  Equation  16  by  taking  into  account  some  spatial  inhomogeneity  of 
the  laser  beam. 

If  the  laser  beam  is  not  uniform,  as  assumed  in  our  derivation,  its  spatial 
energy  density  profile  can  be  represented  again  by  a  gaussian  expression  as 

/  2  An  2  r\2  .  . . 

P  =  PQ  exp  f - g - J  (Al) 

where  is  now  the  integrated  energy  density  at  the  peak  of  both  the  wavelength 

profile  and  the  radial  profile,  r  (  0<  r  <  »)  is  the  radial  coordinate  and  6 

—  —  r 

is  the  spatial  full  width  at  half  maximum  of  the  beam.  In  this  case,  the  inter¬ 
action  of  the  laser  and  the  atoms,  as  given  by  Eg.  (13),  has  to  be  modified 

16  IV 

to  include  the  spatial  dependence  of  the  laser  density.  '  We  therefore 

find  for  the  ratio  of  the  atom  densities  the  following  expression: 


2  An”' 2 


n2  B12  t  V=TlX' 


exp 


f  2 /in  2  ,,  ,  ,V  /  2/in  2  .  V 

p  -  jy —  VVJ  po ,xp  r  -g - V 


"T  h^r) 


2/in  2 


_  - .  -  r  2/ in  2  ,,  ,  41  f  2/ In  l\ 

B'~  T  - expt  “77“  (WJ  po  exp('  6 - /  +  1 


12 


t r  <5A’  **  6A 

which  can  be  written  in  a  simplified  form  as 


n2  B12t  g(5X)pQ  exp  (-cr) 


(A3) 


B12  T  exp  (-  cr)  +  1 


where 


g  (<$A) 


2/C.n  2 

- -  exp 

n  (5x* 


2/ in  2 
fiA* 


'YV2'2 


_  2/in  2 
C  ^  6 

r 


(A2) 
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If  we  now  proceed  as  in  the  text,  introducing  the  saturation  spectral  energy 


density  gi’*en  by  Eq.  (15)  ,  we  obtain 


_/  92  V  exp(-cr)2 _ \ 

°T  \  gi  +  92  ^  exp(-cr)2  +  R  (6X)  / 


(A4) 


where  R  (6X)  = 


P, 


Assuming  the  usual  vertical  monochromator  slit  arrangement  (cylindrical  volume) 


and  a  proportional  variation  of  n^  versus  r  according  to  the  laser  variation, 


the  ratio  in  Eq.  (A4)  can  be  integrated  over  r  as  follows: 


17 


exp(-cr) ' 


2irr  dr 


(A5) 


‘T  J  exp(-cr)  +  R(6X) 
o 


where  the  integration  limit  is  extended  to  infinity  since  the  slit  height  is  much 
larger  than  the  beam  diameter.  By  solving  the  integral,  we  obtain 


t =  ^ £n  * 1  +  Tmr 1 

T  c 


(A6) 


Evaluating  6X  according  to  the  procedure  adopted  in  the  main  text,  we  obtain 
exc 


6X  =  (6)t  +  <<5),~  =  — ^ — -  / 

*  7  nrT  V 


£n 


f-vrrsrj 


(A7) 


which  can  be  written  as 


6X 


6X 


exc 


An  2 


i 


Jin 


(/l+x  -1) 


(A8) 


where  X  5  <PQ/P  )  • 
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f  , 


In  conclusion,  if  we  compare  Eq.  (A8)  with  Eq.  (16),  we  can  see  that  the  spatial 
averaging  effects,  if  present,  tend  to  lessen  the  dependence  of  the  broadening 
of  the  spectral  atom  profile  upon  the  laser  power.  The  effect  will  of  course  be 
more  pronounced  if  the  laser  beam  is  focused  in  the  flame  so  to  reach  high 
irradiances.  However,  even  in  the  presence  of  such  effect,  the  value  of  SX' 
as  obtained  from  Eq.  (A8) ,  is  the  same  as  that  given  by  Eq.  (16) . 


I 
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Table  I.  Experimental  values  of  the  saturation-broadened  half-widths  obtained 
for  Ar-O^-H^  and  the  N2~°2_H2  flam6s-  ^ 


(*) 

All  values  are  within  ±  10%.  Laser  power:  Ca,  5kw;  Sr,  14kw;  Na,  17kW; 
In,  7kW. 
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Table  II. 


Comparison  between  the  theoretical  and  the  experimental 
fluorescence  excitation  profile  half-widths  for  Ar/O  /H 
f 1 ame  s .  ( a ) 


values  of  the 
and  N2/02/H 


Element 

[ ( 6  A  )  / (6A  )  .  ]<b> 

exc  max  exc  min 

Theoretical 

Experimental 

Ar/00/H2 

VV"2 

Ar/02/H2 

V°2/H2 

Ca 

2.5 

2.4 

3.0 

2.2 

Sr 

2.8 

2.8 

2.9 

5.4 

Na 

m 

2.0 

4.3 

1.6 

In 

1.7 

1.4 

2.3 

1.3 

values  are  considered  to  be  within  ±  10%. 

^  (6 \)  refers  to  the  value  obtained  with  the  laser  at  full  power  while  ( 6 A ) 

max  min 

refers  to  that  obtained  when  the  laser  is  attenuated  with  neutral  density 
filters  until  the  fluorescence  signal  is  linearly  related  to  the  laser  irradiance 
(see  values  reported  in  Table  III) . 

(c) 

calculated  according  to  Equation  16  in  the  text. 
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Table  III.  Comparison  between  the  values  of  the  laser  spectral  bandwidth  as 
obtained  by  different  methods 


I 


f 

* 


Element 

1 

1  . 

^  Jb) 

Direct 

Measurement 

Fluorescence  C 
Excitation 

Profile 

(d) 

Saturation 

Broadening 

Ar/02/H2 

VV"2 

*r/o2/H2 

N_/0_/H. 
Z  Z  d 

ca  ; 

| 

0.23 

0.20  j 

0.20 

0.26 

0.21 

Sr 

0.23 

0.42  | 

i 

0.23 

0.46 

0.51 

Na 

0.36 

I 

0.23 

0.40 

0.30 

0.35 

In 

0.24 

0.19 

0.28 

0.26 

values  are  within  +  10%. 


(b) 


values  obtained  by  scanning  the  laser  beam  through  a  1-m  grating  monochromator 

O 

(AX  resolution  =  0.12  A).  Values  are  not  corrected  for  the  instrumental  profile 


values  obtained  by  scanning  the  attenuated  laser  beam  through  the  atomic  vapor 
in  the  flame. 


(d) 


values  calculated  from  the  slope  of  the  plot  obtained  from  Equation  16  in  the 


text. 
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Fluorescence  excitation  profiles  for  sodium  in  the  Ar/O/H,,  flame 
fluorescence  is  shown 


Figure  2D. 


the  resonance 


Fluorescence  excitation  profiles  for  calcium  in  the  flame; 

resonance  fluorescence  at  4227  A;  abscissa:  calcium  concentration 
in  10  vg/mL;  (2)  laser  power:  50W;  6X  =  0.20  A;  (b)  laser  power 

o  6XC 

5  kW;  <5X  _  =  0.44  A. 


I 


Figure  4D.  Fluorescence  excitation  profiles  for  strontium  in  the  ^•*'aine' 

resonance  fluorescence  at  4607  A;  strontium  concentration  in  1  yg/mL; 

(a)  laser  power:  14W;  6Xexc  =  0.23  A;  (b)  laser  power:  1.4  kw; 

6X  =  0.55  A;  (c)  14  kW;  =  1.25  A. 

exc  exc 
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